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MILLIMETER  WAVE  ATTENUATION  IN  MOIST  AIR: 

LABORATORY  MEASUREMENTS  AND  ANALYSIS 

Hans  J.  Liebe 

ABSTRACT 

Experiments  were  performed  with  a  millimeter  wave  resonance-spectrometer 
capable  of  measuring  absolute  attenuation  rates  a  (dB/km)  by  water  vapor  up  to 
saturation  pressures.  Vapor  (e)  and  air  (p)  pressures  were  varied  at  constant 
temperature  and  set  frequency.  Anomalous  absorption  behavior  (i.e.,  high  rates  a, 
extreme  temperature  dependences,  hystereses  in  pressure  and  temperature  cycles) 
could  be  identified  as  being  caused  largely  by  instrumental  condensation  effects. 
Uncorrupted  data  at  138  GHz  display,  in  addition  to  air-broadening  (<*  ep),  a  strong 
self-broadening  component  (<*  e2).  Based  upon  these  results,  a  practical  atmos¬ 
pheric  millimeter  wave  propagation  model  (MPM)  was  formulated  that  predicts 
attenuation,  delay  and  noise  properties  of  moist  air  over  a  frequency  range 
of  1  to  1000  GHz  and  a  height  range  of  0  to  100  km.  The  main  spectroscopic  data 
base  consists  of  48  02  and  30  H20  local  absorption  lines  complemented  by  continuum 
spectra  for  dry  air  and  water  vapor.  Model  input  relies  upon  distributions  of 
meteorological  variables  along  an  anticipated  radio  path.  These  variables  are 
pressure,  temperature,  and  relative  humidity  for  moist  air;  suspended  droplet 
concentration  for  haze,  fog,  cloud  conditions;  and  a  rain  rate.  In  special  cases, 
trace  gas  concentrations  for  ozone  and  carbon  monoxide  and  the  geomagnetic  field 
strength  can  be  added. 

The  research  leading  to  an  MPM  was  reported  in  17  publications  and  is  sum¬ 
marized  in  this  report.  Theoretical  foundations,  experimental  verifications,  and 
typical  radio  path  behaviors  have  been  addressed.  Good  agreement  was  obtained 
between  predicted  and  measured  data  (20  to  337  GHz)  from  both  laboratory  and  field 
experiments.  Implications  of  the  MPM  with  respect  to  older  modeling  codes  and 
needed  future  work  are  discussed  briefly. 


-1- 


1 .  THE  EXCESS  WATER  VAPOR  ABSORPTION  PROBLEM 


Communication  systems  and  sensing  devices  operating  in  the  millimeter  wave 
(MMW)  range  are  attractive  to  military  planners  since  weapon  deployment  has  to  be 
assured  during  periods  of  adverse  weather  and  in  disturbed  environments.  An 
optimum  trade-off  between  atmospheric  obscurations  and  angular  resolution  can  be 
achieved  in  the  frequency  range  between  f  =*  90  and  350  GHz,  where  propagation 
is  restricted  by: 

•  molecular  absorption  due  to  oxygen  and  water  vapor; 

•  mass  concentrations  of  suspended  particles  (haze,  fog,  dust,  etc.); 

•  precipitation  effects  (rain,  wet  snow). 

The  main  problem  studied  during  this  research  effort  was  predicting  correct 
attenuation  rates  by  moist  air, 

a  =  +  av  dB/km 

Dry  air  attenuation  ad,  mainly  due  to  oxygen,  always  makes  a  contribution  which  is 
calculated  with  the  model  detailed  in  [5].  The  line-by-line  summation  in  the 
model  includes  overlap  effects  of  the  60  GHz  absorption  band  and  nonresonant  02 
plus  pressure-induced  N2  absorption. 

Water  vapor  attenuation  ay  depends  on  absolute  humidity  v(g/m3)  which,  in 
turn,  is  limited  by  relative  humidity  RH  (%)  [3].  Its  origin  lies  in  the  absorp¬ 
tion  lines  of  H20  vapor  which  number  in  the  thousands.  In  principle,  ay  can  be 
evaluated  from  a  line-by-line  summation  of  all  line  features,  each  contributing  at 
a  given  frequency  as  prescribed  by  a  line  shape  function  [9].  In  practice,  the 
problem  is  reduced  to  an  evaluation  of  local  line  contributions  a£  to  which  a 
continuum  ac  is  added,  accounting  for  excess  absorption  not  covered  by  the  limited 
line  data  base;  that  is, 

+  ac  dB/km 

The  MMW  frequency  range  has  four  atmospheric  windows  [11]  where  contributions 
by  ac  can  be  as  high  as  90  percent  of  ay.  The  excess  has  been  a  major  source  of 
uncertainty  in  the  correct  modeling  of  moist  air  attenuation  a. 

Basic  physical  principles  underlying  the  attenuation  rate  a  are  difficult  to 
deduce  from  measurements  in  the  actual  atmosphere.  Reliability,  precision,  and 
scale  of  supporting  meteorological  data  compromise  the  quality  of  most  field 
observations.  Controlled  experiments  simulate  atmospheric  conditions,  thus  pro¬ 
viding  test  cases  for  studying  specific  loss  contributions  in  isolation.  Results 
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from  such  laboratory  experiments  led  to  the  formulation  of  a  water  vapor  continuum 
spectrum  ac.  The  continuum  together  with  30  local  H20  absorption  lines  forms  the 
basis  for  an  atmospheric  millimeter  wave  propagation  model.  The  contribution  ac 
has  been  tentatively  identified  to  be  a  far  wing  response  of  the  H20  spectrum 
above  1  THz  [16]. 

2.  MILLIMETER  WAVE  ATTENUATION  BY  MOIST  AIR 

The  most  important  result  of  this  research  is  a  practical  atmospheric  milli¬ 
meter  wave  propagation  model  (MPM)  based  upon  meteorological  input  variables  [14]. 
The  MPM  predicts  attenuation  a  (dB/km)  and  dispersive  delay  B  (ns/km)  up  to  1000  GHz 
Interrelations  between  radio  wave  propagation  and  geophysical  properties  of  water 
vapor  are  a  complex  problem  as  outlined  in  a  comprehensive  overview  [3]  (see 
Appendix).  As  the  work  progressed,  various  theoretical,  analytical  and  experimental 
efforts  followed,  addressing  separate  topics  within  the  framework  set  by  this  review 
Overall,  a  total  of  17  publications  [1]-[17] ,  8  internal  reports  [18] - [25 ] ,  work  by 
the  P.I.  on  4  related  problems  [26]-[29],  and  15  presentations  have  been  produced 
in  order  to  come  to  grips  with  moist  air  as  a  propagation  medium  for  millimeter 
waves. 

2 . 1  Theory 

Theoretical  consultant  for  the  project  was  the  well-known  molecular  physicist. 
Professor  M.  Mizushima  of  the  University  of  Colorado.  He  recalculated  the  rota¬ 
tional  spectrum  of  H20  for  quantum  numbers  J^16  in  the  vibrational  ground  state 
(v  =  0)  and  the  rotational  spectrum  of  02  including  vibrationally  excited  states 
(v  =  0,1,2),  and  applied  the  results  to  predict  atmospheric  transparency  up  to 
frequencies  of  1000  GHz  [9],  [10]  (see  Appendix).  A  detailed  computer  code  is 
available  in  FORTRAN  IV.  Some  new  ideas  with  regard  to  collision-broadened  line 
shapes  were  outlined  [19],  and  contributions  to  excess  absorption  were  estimated 
by  possible,  though  very  weak,  electric  dipole  transitions  in  isotopic  species  of 
N2  and  02  [20]. 

2.2  Modeling 

Modeling  reduces  the  theoretical  complications  of  describing  interactions 
between  MMW's  and  atmospheric  molecules  to  a  parametric  presentation  of  engineering 
utility.  It  runs  the  gamut  from  meteorological  observables  assessing  the  absorber 
population  within  a  radio  path,  through  specific  attenuations,  to  experimental  and 
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theoretical  aspects.  Intricate  "molecule-by-molecule,"  "line-by-line,"  and  "layer- 
by-layer"  calculation  procedures  have  been  converted  into  a  simple  computer  code. 

The  outcome  of  observations  can  be  predicted  provided  the  geophysical  variables 
are  available.  In  general,  agreement  with  experimental  data  at  frequencies  between 
20  and  337  GHz  is  within  limits  of  their  stated  uncertainty  (typically  better  than 
±  10  percent). 

The  following  provides  guidance  to  the  output  of  the  modeling  effort: 
millimeter  wave  propagation  model  MPM,  latest  version  -  [14]; 

program  descriptions  of  the  MPM  -  [18],  [14],  [12.’’ 

local  line  model  (30  H20  and  48  02  lines)  -  [4],  [E 

water  vapor  continuum  absorption  (see  Appendix  A4)  -  [11],  [13],  [15],  [If 
simplifications  of  the  MPM  -  [12],  [1. 

applications  of  the  MPM  -  [7],  [8],  [11],  [13],  [14],  [23],  [24],  [27]-[29j. 

A  typical  graphical  output  of  the  MPM  is  shown  in  Figure  1.  When  examining  the 
attenuation  rates  a,  as  listed  in  Table  1  for  "window"  and  "line"  (185  GHz)  fre¬ 
quencies,  on  a  per  g/m3  concentration  (v)  basis,  one  notices  in  Figure  2  striking 
differences  between  the  old,  Gaut-Reifenstein-based  model  [11]  and  MPM  predictions. 
Discrepancies  at  warmer  temperatures  T  are  caused  mainly  by  a  v2-dependence  of  the 
continuum  ac  (see  Eq.  29  -  p.  A27)  which  are  amplified  by  the  strong  temperature 
(  *  T17)  dependence  of  vs(100%  RH). 

2.3  Experiment 

The  experiment  consists  of  a  millimeter  wave  resonance  spectrometer  and  an 
integrated  humidity  simulator.  The  instrument  is  capable  of  measuring  absolute 
electromagnetic  properties  (attenuation,  dispersion,  refraction)  of  moist  air  for 
relative  humidities  between  0  and  100%  RH  [16]  (Appendix).  An  insulated  box  con¬ 
tains  a  high-vacuum  stainless  steel  vessel,  which  houses  a  temperature-controlled 
mini-lake  (10  cm  across)  and  the  radio  test  path.  MMW  signals  around  35  and  140  GHz 
have  been  employed  over  equivalent  free-space  path  lengths  between  350  and  150  m; 
temperatures  are  controlled  to  better  than  1/100  of  a  degree  Celsius;  pressure 
ranges  over  seven  orders  of  magnitude  (103  to  10"4  mb);  and  relative  humidity  is 
varied  between  0  and  100  percent.  Data  acquisition  is  computerized,  with  a  choice 
of  operational  modes  covering  time  scales  from  10  seconds  to  many  days.  The  reson¬ 
ance  signal  that  carries  the  absorption  information  is  detected  with  a  digital 
waveform  analyzer. 

Mechanical  design  details  are  given  in  [21]  and  results  of  tests  on  tem- 
perature-vs-humidity  performance  are  presented  in  [22].  The  millimeter  wave 
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Figure  1.  MPM  predictions  of  specific  attenuation  a  (dB/km)  for  dry  and  moist 
air,  including  fog  (about  300  m  visibility)  and  rain  events  [14]. 


sources  [that  is;  klystrons  centered  at  35  GHz  (OKI  35V11),  at  73  GHz  (OKI  70V11A), 
and  at  139  GHz(VARIAN  VC  71 5S),  and  a  Gunn  oscillator  centered  at  92  GHz  (MYCRON)] 
operated  CW  in  a  phase-locked  mode  with  10  kHz  resolution,  or  frequency-stepped 
under  program  control.  In  addition,  a  free  running  swept- frequency  mode  was 
provided.  A  tripler  (MYCRON,  6P2  diode)  generated  73  x  3  *  220  GHz  power  in 
amounts  between  0.5  and  2  mW.  Performance  curves  are  given  in  Figure  3.  The 
tripler  varactor  diode  failed  repeatedly  (6  times,  so  far).  Consequently,  we  have 
not  yet  been  able  to  perform  moist  air  attenuation  measurements  at  220  GHz.  Some 
circuit  details  of  the  phase-lock  loops  are  shown  in  Figure  4. 
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Table  1. 

Predicted  specific  attenuation  a  (dB/km)  of  sea  level  air  (p  =  1013  mb) 
at  various  relative  humidities  RH  (%)  and  temperatures  T(K)  for  five 
millimeter  wave  frequencies  f  (GHz)  [14]. 


Specific  Attenuation  a 

f 

...  T 

1 00%  RH 

75%  RH 

50%  RH 

25%  RH 

0%  RH 

35  GHz 

K  g/m3 

■mi 

dB/km 

300  25.49 

mgmMm 

0.275 

0.172 

0.090 

0.028 

290  14.31 

0.153 

0.105 

0.031 

280  7.65 

0.119 

0.095 

0.052 

0.034 

0.079 

0.068 

0.057 

0.047 

0.038 

95  GHz 

300 

2.239 

1.475 

0.853 
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Figure  3.  Power  curves  of  the  220  GHz  synthesizer:  (a)  stepped-frequency  mode  under  program  control  through 
resonances  of  two  wavemeters  centered  at  217.95  and  222.2  GHz;  (b)  reduced  resonance  curves;  and 
(c)  retraces  of  the  resonance  center  for  210  (=  3.5  hrs)  and  75  (=  1.25  hrs)  scans  each  in  50  steps. 


Parametric  studies  of  the  attenuation  rate  a  were  performed.  Experimental 

variables  are  frequency  f,  total  pressure  P,  temperature  T,  and  humidity  v(RH). 

Absolute  attenuation  due  to  water  vapor  and  moist  gas  mixtures  has  been  measured  up 

to  saturation  pressures  (RH  =  100%)  and  void  of  the  instrumental  (condensation) 

errors  plaguing  earlier  measurements.  Hysteresis  effects  in  cycles  of  relative 

humidity  (Fig.  11,  p.  A13)  indicated  condensation  on  the  spectrometer  mirrors  and 

★ 

resulted  in  highly  erroneous  attenuation  values.  A  large  body  of  earlier  data  had 
to  be  disqualified  because  of  this.  Slight,  localized  heating  (+  1°C)  of  the 
mirrors  and  a  double-sealed  coupling  iris  cured  the  problem. 

Laboratory  results  in  support  of  the  MPM  are  reported  in  [1],  [6],  [13], 

[15],  [16],  [17],  [21],  [22],  [24],  and  [26].  Measurements  at  35  GHz  [26]  proved 
to  be  inconclusive  since  moist  air  attenuation  a  is  relatively  weak  (Table  1). 
Measurements  at  96  GHz  were  performed  with  a  well  calibrated  horizontal  radio  path, 

27  km  long  and  operated  in  Boulder,  CO  (h  =  1.5  km)  [14],  [24d].  Most  of  the  data 
in  support  of  the  new  continuum  ac  were  taken  at  138  GHz  [16]. 

2.4  Needed  Future  Work 

Two  shortcomings  of  the  propagation  model  MPM  lie  in 

•  the  missing  confirmation  of  the  physical  origin  for  the  water 

continuum  <*  and 
c 

•  an  oversimplified  treatment  of  opaque  atmospheric  conditions 
caused  by  haze,  fog,  clouds,  dust,  smoke,  etc. 

Additional  parametric  studies  of  the  water  vapor  continuum  absorption  for  self- 
and  air-broadening  should  be  performed  at  high  relative  humidities  (RH  *  80  to  100%), 
different  temperatures  (T  =  270  to  320  K),  and  total  pressures  (P  =  100  to  1500  mb) 
at  the  two  window  frequencies,  f  =  140  and  220  GHz  (see  Table  1).  Artificial  aerosol 
populations  of  known  chemical  composition  and  concentration  can  be  added  to  study 
their  growth/evaporation  behavior  as  a  function  of  humidity. 

The  pressure-broadened  water  vapor  absorption  line  centered  at  183.3  GHz 
might  be  included  in  the  parametric  laboratory  studies  to  obtain  accurate  strength, 
width,  shift,  shape  data,  which  are  particularly  important  to  remote  sensing 
appl ications. 

★ 

We  suspect  that  this  source  of  instrumental  error  also  affected  results  reported 
by  other  workers,  thus  contributing  to  the  many  inconsistencies  of  reported  water 
vapor  excess  absorption  [3]. 
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Proceeding*  oC  the  International  Workshop 
on  Atmospheric  Water  Vapor 
Held  in  Vail  Colorado. 

September  11-13,  1 97* 


PREFACE 


Thu  volume  contain  the  technical  procccdtnp  o(  the  International  Workshop 
on  Atmoaphcnc  Water  Vapor  held  in  Vail.  Colorado.  September  11-13,  1979. 

The  workshop  wat  orpnued  to  bring  together  mcarchera  in  various  relaud 
fields  to  aiMSi  the  present  suie-of-the-art  m  measurement.  modeling,  and 
Application  of  atmospheric  water  vapor  properties  and  to  identify  those 
importaaa  problems  in  which  further  effort,  in  the  near  future,  it  needed  for 
better  understanding  of  the  atmosphere  itself,  and  of  electromagnetic  propa¬ 
gation  through  the  atmosphere.  One  hundred  and  eleven  scientists  from 
universities,  government  agencies,  research  laboratories,  and  industry  in  the 
United  States.  England.  France.  Germany,  and  India  attended  the  Workshop  in 
which  fatly  five  pipers  were  prevented  Complete  texts  of  thirty  nine  of  these 
papers  and  tome  selected  discussions  are  included  in  this  volume. 

Atmospheric  water  vapor  plays  an  important  role  in  many  atmospheric 
processes,  such  at  those  related  to  weather,  climate,  atmospheric  chemistry, 
radiative  transfer,  dectrooptie  imaging,  and  remote  te rising  Knowledge  of  the 
global  distnboltoo  of  water  vapor  and  its  transport,  optical  properties 
(absorption  and  emission)  in  the  infrared  and  microwave  regions,  and  influence 
on  electromagnetic  propagation,  are  of  particular  interest  to  a  wide  range  of 
rciemific  disciplines. 

The  workshop  program  was  divided  into  seven  piper  presentation  sessions, 
one  open  discussion  and  one  panel  session.  Thirty  invited  and  fifieea  con¬ 
tributed  papers  covering  the  various  a i pacts  of  atmospheric  water  vapor  were 
presented  to  the  seven  sessions  under  the  excellent  guidance  of  the  union 
chairmen-  as  indicated:  S.  A.  Cough.  Optics  and  Spectroscopy  of  Water  Vapor. 
H.  J.  Lube,  Microwave  and  Millimeter  Wave  Phenomena:  G.  M.  Reach. 
Geoastrophysicai  Applications.  P  M  Kuhn.  7n  55 u  Measurements  of  Water 
Vapor.  V.  £  Derr,  Remote  Sensing  of  Water  Vapor  F.  Baer.  Meteorology  of 
Water  Vapor.  V.  Mohnen.  Atmospheric  Chemistry  and  Microphytics  of  Water 
Vapor.  The  papers  discussed  the  current  state  of  knowledge,  as  well  atthc  results 
of  the  latest  investigations  in  specific  areas  of  research.  The  moderators  for  the 
open  and  panel  discussion  sessions  were  A.  Gebbie  end  V.  Mohnen.  respectively. 
Ample  time  was  allowed  For  discussions  following  each  paper  and  in  the 
dtscaseion  sessions.  Discussions  wen  recorded  and  the  transcripts  post-edited. 

To  ensure  proper  representation  of  major  disciplines  involved,  a  Technical 
Program  Comnwice.  composed  of  the  Following  individuals,  was  set  up-  A. 


Deepak.  Institute  For  Atmospheric  Optics  and  Remote  Seining  (FF  AORS). 
Associate  Chairman;  L.  H.  Ruhnltr.  Natal  Research  Laboratory  (N’RL), 
Cochairman;  T.  D.  Wilkcrson.  University  o(  Maryland.  Cochairman;  K. 
Jayaweera.  National  Science  Foundation  (NSF);P.  M.  Kuhn.  National  Oceanic 
and  Atmospheric  Administration  (NOAA);  G.  M.  Reach.  Jet  Propulsion 
Laboratory  (JPL);  J.  Rogers.  Federal  Aviation  Administration  (FAAk  and  S 
Titford.  National  Aeronautics  and  Spare  Administration  (NASA)  Hesd- 
quaners. 

The  workahop  was  cosponsored  by  Naval  Research  Laboratory.  Naval  Air 
Systems  Command.  Office  oF  Naval  Research  (ONR).  National  Aeronautics 
and  Space  Administration.  National  Oceanic  and  Atmospheric  Administration. 
Federal  Aviation  Administration.  National  Science  Foundation.  University  oF 
Maryland,  and  the  Institute  For  Atmospheric  Optics  and  Remote  Sensing. 

The  editors  wish  to  acknowledge  the  enthusiastic  support  and  cooperation  of 
the  members  of  the  Technical  Program  Committee,  session  chairmen,  speakers, 
and  participants  foe  making  this  a  very  mmulalmg  and  highly  beneficial 
workshop  For  everyone.  Special  (hanks  arc  due  the  authors  for  their  cooperation 
in  enabling  a  prompt  publication  of  the  workshop  proceedings.  It  is  s  pleasure  to 
acknowledge  the  valuable  issistsncc  of  IFAORS  staff,  especially.  Sherry  Allen. 
Sue  Crofts,  and  June  Ewing,  in  the  smooth  organization  of  the  workshop,  and 
H.  Mulcahy  and  M.  Goodwin  in  preparing  and  typing  of  the  manuscripts. 

The  editors  hope  this  volume  will  be  a  useful  contribution  to  the  field  of 
atmospheric  water  vapor  and  serve  ts  a  reference  volume  for  many  years  to 
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Atmospheric  Water  Vapor. 

A  Nemesis  for  Millimeter  Wave  Propagation  1 
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Oxygen  absorption  .... 
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Millimeter  vsvts  offmz  an  attractive  way  of  salving  uni  qua 
system  problem  because  of  their  ability  to  penetrate  clouds, 
smog,  dust  at  fog.  This  make*  them  s  logical  choice  over  electro- 
optical  devices  for  adverse  weather  applications .  Spectral  lines 
of  oxygen  and  water  vapor  ultimately  limit  the  atmospheric  trans¬ 
parency/  hence ,  mat  applications  operate  between  the  absorption 
lines  in  four  window  regions  < that  is,  24  to  4$,  72  to  110*  129 
to  140,  and  200  to  260  CM*) .  Observations  have  established  the 
esistence  of  excess  water  vapor  absorption  (EWA),  which  is  most 
evident  xn  these  windows.  Excess  implies  that  the  effect  is  not 
related  to  the  known  spectral  properties  of  the  water  molecule. 

EWA  is  found  to  increase  in  nontrivial  manner  with  humidity  and 
the  discrepancies  can  be  as  large  as  a  factor  of  10.  Several 
groups  (most  notably  at  the  Appleton  Laboratory.  UK) ,  have 
gathered  evidence  of  EM*  from  laboratory  and  field  observations 
and  brought  forward  hypotheses  to  account  for  the  date. 

Qualitative  explanations  are  based  on  the  assumption  that  watsr 
molecules  in  moist  air  form  larger  molecules  with  a  dimer  being 
the  first  step  in  a  series  of  stable  species.  Hydrogen  bonding, 
ion  attraction  and  attachment  of  the  polar  MOO  to  foreign 
par  tic  lee  (aeroeol  growth )  arm  the  ordering  forces  considered  in 
the  clustering  process.  An  assessment  of  the  current  EWA  picture 
will  be  given  and  avenues  of  research  attacks  are  discussed  to 
solve  the  enigma  in  the  quantitative  description  of  the  inter¬ 
action  between  millimeter  waves  and  moist  air. 


1Work  was  partially  supported  by  the  O.S .  Any 
under  apo  JO-79. 
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1 .  INTRODUCTION 


Atmospheric  water  in  both  vapor  and  liquid  state*  1*  the 
maior  deterrent  to  an  unrestricted  exploitation  of  propagation 
of  millimeter  and,  sore  so,  of  infrared  wavelength*.  Tor  moat 
applications,  the  operation  of  ground-based  systems  is  limited 
to  seven  wlndw  regions  wl  to  W7,  these  being  the  gaps  between 
molecular  absorption  lines  and  bands < 


Absorption  feature 


Region  window  Range 
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the  sain  attraction  of  millimeter  wave  systems  is  their 
ability  to  penetrate  the  somewhat  opaque  atmosphere  (base,  fog, 
elouds,  dust,  smoke,  light  rain)  undsr  circumstances  in  which 
electro-optieal  and  lnfrsrsd  systems  normally  fall.  Accurate  and 
detailed  knowledge  of  atmospheric  transmission  la  esaantlai  to  an 
evaluation  of  the  advantages  of  millimeter  waves  over  the  shorter 
wavelengths. 


One  objective  of  this  paper  is  to  focus  attention  on  unsolved 
problems  in  the  construction  of  a  valid  model  for  atmospheric 
millimeter  wave  transmission,  subject  to  the  following  conditional 
Frequency,  v  •  10  to  1000  GHx.  with  special  emphasis  on  the  CHF 
range,  v  •  30  to  300  GHtt  altitude,  h  •  0  to  30  lorn  and  relative 
humidity,  RH  •  0  to  <,  1009.  The  radio  engineer  usee  the  optical 
term  'clear*  for  moist,  but  un saturated  air  (RH  <  100%) ,  and 
ignores  hare  conditions.  Although  many  constituents  contribute 
to  the  total  atmospheric  attenuation  rate*  a  MB/ km) ,  absorption 
in  the  window  is  dominated  by  water  vapor  and  is  of  greatest  con¬ 
cern  for  practical  situations.  At  present,  describing  water 
vapor  absorption  mechanimas  and  rslating  than  to  measurable 
quantities  is  partly  an  empirical  matter  and  lacke  credibility  in 
a  general  sense. 


The  following  interaction  effects  between  millimeter  waves 
end  air  are  Mil  documented! 


110 


(9)  Complete  line  parameter 
pi la t ion  for  all  etas- 
spheric  molecules  .  . 
known  as  the  ATGL  Tape* 

H1TRAN  Model  and  rASCOOE  .  10.  14.  25,  29 

(h)  Scattering  loss  by  particles  33*  32.  3S-37 

(i)  Dielectric  loss  of  watsr  .  .  25,27,  13,  102-106 
Other  physical  effects,  which  have  been  taken  into  account  in 
infrared  transmission  models  (9),  (20),  and  (23)  are: 

(k)  Submicron  aerosol  particles  .  38,  41.  42 

(l)  Growth  of  subeiicron 

particles  in  moist  air  .  .  34-37 
Of  a  more  speculative  nature  are  considerations  of 

spectra  .  6*  fl,  21.  43*  58.  64.  66, 
69.  61,  86,  89.  109 
,  n>2  .  .  21,  22.  44 

(o)  Water  ion  activity . 28.  45 


The  water  polymer  hypotheses  have  the  inherent  weakness  that 
except  for  the  dimer  (43).  neither  definite  structural  observa¬ 
tions  have  been  made  nor  millimeter  wave  spectre  calculated.  Zn 
the  cases  (m)  to  <o) .  production  and  growth  mechanisms  (for 
example,  ion-induced)  as  well  as  size  and  nimber  density  distri¬ 
butions  have  not  been  Identified  for  the  natural  atmospheric 


Measured  absorption  in  the  windows  W1  to  W7  could  provide 
clues  for  ths  understanding  of  which  basic  physical  mechanisms 
are  missing  In  current  models,  new  instrumentation,  such  as  the 
dispersion  spectrometer  (63.65)  and  the  saturation  hygrometer 
(107) ,  might  prove  helpful  in  future  investigations.  Available 
experimental  results  on  absolute  attenuation  rates,  path  trans¬ 
mittance*,  and  thermal  emission  from  laboratory  (46-71)  and  field 
(73-171)  observations  in  cloudless  moist  sir  arc  generally  higher 
than  values  predicted  on  the  basis  of  molecular  absorption  alone. 
It  is  a  fact,  even  after  allowing  for  difficulties  in  measuring 
the  highly  variable  water  content  quantitatively,  that  the 
absorption  by  water  vapor  in  the  atmospheric  transmission  windows 
wi  to  W?  iB  not  completely  understood. 


Section  II  of  thie  paper  presents  the  framework  for  s  radio 
path  modeling  scheme  of  which  details,  such  as  measurements  and 
validation,  have  been  published  (15.17,19).  The  model  serves  as 
a  basis  for  the  definition  and  discussion,  in  Section  XV,  of  the 
problem  of  excess  water  vapor  absorption  (EWA)  in  the  light  of 
experimental  data.  Section  III  is  devoted  to  a  brief  presenta¬ 
tion  of  relevant  physical  properties  of  atmospheric  water  vepor. 


II.  EHF  RADIO  PATH  MODELING 


Millimeter  waves  traveling  through  the  cloudless  (RH  _<  100%) 
atmosphere  suffer  both  deterministic  and  random  variations  in 
amplitude  and  phase.  Radiation  is  absorbed  and  refracted  by 
gases  and  submicron  particles.  This  lntsraction  is  modeled  with 
the  intent  to  predict  the  foliating  frequency-dependent  propaga¬ 
tion  effects:  (a)  absorptive  loss  of  coherent  radiation. 

(b)  time  of  propagation  between  two  points,  (cl  rsfr active  ray- 
bending  and  ducting,  (d)  generation  of  incoherent  noise,  and 
(e)  scintillations  due  to  random  fluctuations  of  the  mcdii*  in 
space  and  time.  The  array  of  mathematical  and  empirical  expres¬ 
sions  used  to  compute  these  various  effects  is  called  the 
Propagation  Model.  Such  a  parameterization  scheme  depends  fore¬ 
most  on  spectroscopic  information  about  kind  and  number  density 
of  absorber  species  and  on  their  distribution  within  the  path 
volume. 


A.  The  Propagation  Model 


Complex  refract ivity  N(ppm)  is  a  convenient  macroscopic 
measure  of  the  Interaction  between  millimeter  wave  radiation  and 
the  individual  absorbers  In  moist  air.  A  value  of  N  accounts  for 
the  effectiveness  and  number  density  of  a  particular  absorber 
population.  Later  these  dependences  are  formulated  in  terms  of 
measurable  quantities.  At  this  point,  N  is  assumed  to  be  known 


-A2- 


* , *  ■  ■  %*  *r •  /■  .  -  »*•  a*'  » ■-  v 


^axi  '  •  i  **■*  ‘  1  I 


A  NEMESIS  F0«  MILLIMETER  WAVE  PROPAGATION 


S 


Ir';’’ 


and  the  basic  relations  that  determine  the  electromagnet ic 
behavior  of  a  radio  path  ara  ««tt*n  down.  Tha  comple*  refrac¬ 
tivity  of  a  gaseous  medium  in  pp*’' 


D(V)  *  )  N*|V) 


(1) 


consist a  of  thraa  components;  namely*  tha  frequency  independent 
refractivity  Kq  plus  various  spectra  of  refractiva  dl spars  ton 
D(v)  and  absorption  N“(ul.  Tha  atmosphere  is  charactarisad  as  a 
linear  network  with  both  passiva  and  active  properties.  The 
amplitude  and  phase  response  of  a  plane  radio  wave  travel  Inf  tha 
distance  L(k»)  and  having  an  initial  field  strength  Eq  is 
described  by 


t  •  E  exp  (ft.) 


(3) 


10  > 


(3) 


r  •  j(2»v/e)U  ♦  n 

is  the  propagation  constant  of  the  intervening  medlue,  c  bainf 
tha  spaad  of  lifht.  Usually,  raal  and  imaginary  parts  of  N  are 
separated  and  expressed  as  the  power  attenuation  rata  in  dB/Tos 


a  -  o.lflJOvDMN) 

and  the  phase  daisy  rate  in  radian/km 
*  -  0. 0209vRE (N) 

Tha  frequency  U  Is  in  gigahertz  iZVit)  throughout  the  papot. 


(4) 


The  attenuation  rate  a  is  the  more  familiar  quantity  in 
atmospheric  wave  propagation.  The  phase  rate  e  must  be  consid¬ 
ered  when  it  varies  spatially  (for  example,  radar  pointing 
accuracy,  long-baseline  interferometry,  maxima*  dish  size  for 
reflector  antennas,  etc.)  or  with  frequency  (for  example,  band¬ 
width  limitations  of  a  cewwunication  channel).  Characteristics 
of  short,  horizontal  radio  paths  may  be  approximated  by  an 
average  value  of  N.  more  general  path  geometries,  such  as  a 
ground - to- sa tell ite  link,  are  treated  by  dividing  the  path  into 
segments  having  quasi -constant  H- values  and  suawing.  The  omsmi- 
1. stive  behavior  between  the  ray  points  si  and  sj  is  expressed  in 
dB  by  the  total  attenuation 


!• 


or  by  the  transeittance  (multiply  by  1/10  log  e  •  0.23026  to 
convert  dl  to  Np 


t  •  exp (-0.230*) 


(7) 


The  medium  becomes  transparent  when  T  approaches  one  and,  on  the 
other  hand,  opaque  for  t  *  0. 


The  total  phase  change  for  the  same  path  is  given  in  radians 


*fj 

y(s)  ds 


<•) 


which  translates  into  the  travel-time  in  ns  of  the  wave  by 

7  x  6/2W  (9) 


The  path  differential  ds  is,  in  practice,  an  lncrtswnt  ds  over 
which  n  is  quasi- constant*  and  depends  upon  the  altitude,  hi  the 
starting  angle,  9,  from  the  zenith  in  the  case  of  a  slant  pathi 
and  refractive  bending  of  the  ray  (Snell's  law)  due  to  gradients 
de/5s  (19). 


Absorption  by  the  atmosphere  causes  emission  spectra.  Each 
unit  of  volume  maintains  thermal  equilibrium  with  its  environment 
via  collisions;  hence,  the  path  element  ds  radiates  an  equivalent 
blackbody  emission  T(s)  a(s)  ds,  which  is  reduced  by  the  trans¬ 
mittance  Tls)  along  the  path  of  observation.  The  resulting 
brightness  tc^>ersture  in  K 


T(s)  a<a)  T (s)  ds 


(10) 


is  either  less  than  or  egual  to  the  ambient  temperature  T. 
Equations  (6),  (9),  and  (10)  constitute  the  key  by  which  perfor¬ 
mance  limitations  of  EHT  systems  operating  over  clear-air  propa¬ 
gation  paths  may  be  evaluated.  A  transfer  function  exhibiting 
constant  amplltua  (A  -  constant),  frequency- linear  phase  delay 


(Dfv)  -  0),  and  no  noise  (Ts  -  0)  implies  ideal  channel  behavior. 


A  broadband  signal  occupying  a  frequency  interval  dv  is  distorted 
by  the  deterministic  spectra  of  D(v)  and  N*(v|>  in  addition,  TB 
imposes  detection  limitations.  As  a  bonus,  the  emission  spectra 


Tg(v)  afford  opportunities  to  sense  remotely  the  state  of  the 


atsvsphere  by  peesive  radio 
40). 


trie  means  (for  example,  Ref.  39 


B.  The  EHF  Refractivity  of  Moist  Air 


The  physical  state  of  moist  sir  Is.  described  by 


dry  sir  prsssurs  (1  kPa  •  10  mb)  P  ,  kPs 

rslativs  invsrss  temper at ue  <T  m  ri  t  •  300/T 
vatsr  vapor  partial  pressure  •  .  kPa 


Calculation  of  the  frequency-ind»q  indent  refractivity  in  ppm  in 
Eq.  (1)  is  straight forward  (47-50 » 


2.509  p  t  ♦  (41.6  t  ♦  2.39)  et 


(12) 


watee  vapor  refractivity  la  about  16  times  more  effective,  on  a 
per  molecule  basis,  than  dry  air  m  generating  propagation 
phenomena  such  as  daisy,  ray  bending,  ducting,  scintillations, 
etc. 


The  dispersion  contribution  in  ppm 
2 


and  the  absorption  spectrin  m  ppn 
W(v)  -  l.  (a  F") .  ♦  M“  ♦  n- 


rsquirs  further  elaboration, 
spectra  are  of  two  types: 


Frequency -dependent  molecular 


1.  Line  spectra  of  absorption  SF*  and  of  dispersive  refrac¬ 
tion  sr' ,  having  strength  S  in  units  of  kHz  and  shape  factors  r* 
and  F”  in  units  of  GHz*1;  the  suns  over  i  consider  millimeter 
wevs  lines  (see  Tables  1  and  2)  of  Oj  (1  •  2  to  45)*  and  HjO 

(1  *  46  to  74).  Spectra  of  the  trace  gases  O3 ,  CO,  M^,  SOj, 
mi),  etc.,  are  neglected  (10,17.29,31.110). 

2.  Continuum  water  vapor  spectnsn  due  to  far-wing 
contributions  of  very  strong  infrared  lines.  A  third  term,  N*, 
which  is  not  fully  understood,  was  sdded  to  the  absorption  to 
account  for  contributions  other  then  those  of  the  rotational 
watar  vapor  line  epeetri—  (so#  Section  IV) . 


Coe— on  to  each  spectroscopic  feature  is  an  intenslty-agalnst- 
frequency  distribution  function,  tha  shapes  F* (v)  and  F~(v). 

For  h  <  20  km,  the  shape  functions  are  (17)a(f/) 


(v  -v>  ♦  yt 
—  o 

L - “ - 


(v  *  V)  ♦  Yl 


(IS) 


The  nonraaonant  oxygen  spectrum, 
Rat.  17. 


i  ■  J  js  discussed  in 


FIGURE  2.  Attenuation  raee  z  at  sea  levei  (h  •  0  km.  e  * 

1.015  kfa)  over  the  frequency  range,  v  •  30  CM*  to  JO  r*t  display¬ 
ing  the  anvalopa  for  masimum/ninimum  values  of  the  rotational 
voter  aapor  spectrum.  Calculate  by  using  the  AFGL  Tepe  (20)  . 
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r" 


Y  -  (v  -  v)  I 
o 

iv  -  v>*  ♦  y2 

o 


-  (V  t  V)  I 


(16) 


Introducing  the  spectroscopic  ps rameters  of  center  frequency  v0, 
width  Y.  end  overlap  interference  I.  F*  and  F*  are  in  units  of 
l/GH*. 

For  h  >  20  Jus.  only  isolated  resonant  lines  are  present, 
spread  over  a  megahertz  frequency  scale.  Equation*  (15)  and  (16)  Ft 
reoUieJ  Lorentzian  shapes 

r;  .  (v.  -  vl/Hv,  -  v)J  ♦  yJ)  ( ID 

L  O  O 


F£  -  Y/((Vfl  -  V)2  ♦  Y2) 


(18) 


Very  accurately  Measured  trfrsctlvity  can  serve  as  a  test  for  th*» 
low  frequency  response  of  theoretical  line  shapes  applied  to  the 
rotational  HjO  spectrum  if  one  is  sure  of  the  spectroscopic  data 
baa*.  A  more  detailed  discussion  of  Eqs.  (Jf)snd  (27)  follows  In 
Section  IV. 


4.  Tempera Cure  Behavior.  A  closer  look  at  the  spectral  line 
resnnnse  elds  In  underatending  the  temr«r*ture  dependence  of  EHF 
attenuation  end  pheee  dispersion  rates  in  moist  air.  Lin*  lnten- 
sitlea  are  senaitlv*  to  temperature  variations  at  various  rates. 
The  peek  line  intensities  in  ppm  (Eq.  (19))  are  described  by 


_! _ S£ _ 

°2  U^bjJlp/lp  ♦  !•!•»)  t2,1e*pU2U  -  t>) 

BjO  (b1/bJ)ls/(p  ♦  4. Be)  J  t2*%*p(bj  (1  -  t)j 


(29) 


where  F£  and  F£  are  in  unite  of  1/GHz.  Peak  dispersion  at 
v  -  v0  »  Y  and  maximum  absorption  at  v  •  v0  are  given  in  ppm  by 

10.  - s/* 

H*  -  SA  d»> 

For  h  >  40  km.  a  further  decrease  in  pressure  converts  the 
pressure-broadened  Lorentzian  into  a  Doppler- broadened  Gaussian 
ahape  with  a  different  width  in  kHz  of 

y  -  «.jo  v  /«5t  i JO) 

D  O 


The  line  intensities  ere  independent  of  tarpereture  when 


t  -  1.1  1.3  1.5  2 

°2 

#2  -  2.11a  t/(t  -  1) 

2.00  1.97  1.70  1.46 

(see  Table  1) 

*2° 

b2  •  2. 9fn  t/(t  -  1) 

2.76  2.54  2. 35  2.02 

(see  Table  2  ) 

Lines  with  values  lower  then  those  prescribed  by  Eq.  (30) 
increase  in  intensity  whsn  the  temperature  drops  (t  y  1)  and  vice 
versa.  A  useful  approximation  at  a  fixed  frequency  is  given  by 


where  a  is  the  molecular  weight  (e.g..  Oj,  m  -  32.  t  •  1. 

V0  -  60  GRZt  YD  ■  65.8  kHz). 

The  convolution  of  Lorentzian  and  Gaussian  shape  functions  is 
called  the  Voigt  profile,  which  is  governed  by  the  parameter 
Y/Yd-  The  Voigt  profile  is  appropriate  when  this  ratio  Celle  in 
the  range  between  *  10  and  0.1.  Numerical  evaluation  of  the 
height-dependent  complex  Voigt  function  requires  considerable 
cosgmtational  effort  (19). 

1 .  Micro wwve  Spectrum  of  Oxygen  (Oj-HSJ .  The  Dj  molecule 
has  electronic,  vibrational,  and  rotational  energy  levels  with 
transitions  causing  spectral  lines  from  the  optical  to  the  EHF 
range.  The  EHF  lines  are  fine  structure  transitions  between 
rotational  triplet  energy  levels.  All  pertinent  information  on 
the  02-W  is  tabulated »  that  is.  the  line  parameters  have  been 
reduced  to  five  coefficients  a}  to  a^.  which  are  listed  in 
Table  X  together  with  the  center  frequencies  v0  (17.19). 
Dependence  on  atmospheric  properties  is  sxprsssed  by 

1  ■  I,  p  t!  npltjll  -  til  (211 

Y  •  *j  (P  *  1.1.1  t°'9  (11) 

l  »  «4  p  t**  (21) 

where  S  and  Y  are  in  kHz  and  GHz.  respectively.  Hater  vapor 
influences  the  02-ms  through  Bq.  (22) . 

2.  Microwav  Spectrum  of  Hater  Vapor.  From  the  ATCX  line 
parameter  compilation  (10) ,  1838  HjO  lines  wars  reed  up  to  31  THx 
when  £  -uxisa  intensity  (v  •  v0)  cutoff  of  2  ■  10“ 3  <31/ km 

(3QQ  K)  is  applied.  FOr  the  EHF  range.  29  line*  (see  Table  2) 
must  be  considered  explicitly  end  the  remaining  ones  are  lushed 
into  e  continuum  fax-wing  contribution.  Mi  individual  line  is 
described  by 


bj  •  t5'5  .xplbjll  -  t)) 

(24) 

bj(p  ♦  «.ao  *>t°'s 

(25) 

-  o 

where  Sy  end  Yv  ere  in  kHz  end  GHz.  respectively. 

3.  Continuum  Hater  Vapor  Spectrum.  The  remaining  1809 
(that  is.  1838  -  29)  HjO  lines  make  far-wlng  contributions t  they 
er*  fitted  by  (17) 

11“  ^  *  1.9  x  10“*  pet4  (V/30)  (26) 


a(t>  -  at 300  k>  ty  131) 

where  the  exponent  y  is  obtained  by  fitting  nodal  date. 

C.  The  Radio-Path  Model 

The  path  la  asoimmd  to  be  in  a  spherically  stratified  atmo¬ 
sphere  In  which  each  layer  is  homoaeneoos.  The  integral 
expressions  (Eqs.  (6).  (8)  end  (10))  are  evaluated  by  the  isyer- 
hy-layer  method  (19,29)  by  using  numerical  integration  techniques 
since  closed-form  solutions  ere  lacking  because  of  the  complicated 
temperature  height  profile. 

Mean  conditions  are  modeled  by  the  0.5.  Standard  Atmosphere 
1976,  and  in  altu  date  from  radiosonde  ascents  can  be  pr  oararmed 
in  directly.  Whatever  description  of  the  atmosphere  Is  employed, 
model  or  date.  It  is  converted  in  the  computer  into  n- layers 
each  having  an  assigned  set  of  values  of  p(h),  e(h)  <«,,  t(h). 
^luea  for  a  are  calculated  by  usinq  Eqs.  (41,  (14),  (161  and 
(21)  to  (25)  and  form  by  applying  Eqa.  (5),  (12),  (13.)  (15)  and 
(21)  to  (25),  by  the  standard  Jinc-hy-Jine  superposition  (29)  and 
adding  the  continuum  toq,  (27). 

The  radio  path  is  assumed  to  be  a  ray  (that  is,  plane  wsve 
case),  and  the  general  procedure  is  to  calculate  the  local  N(h) 

(1)  for  the  programmed  altitude  grid  reporting  p,  e,  t  and  to 
store  n0,  D(V),  and  w"fv)  separately  on  file.  The  ray  starts 
at  a  surfacs  height  with  *n  elevation  angle  8  and  is  guided 
through  the  inhomogeneous  medium  until  it  reaches  the  Intended 
final  height  h}.  Refraction  Is  most  pronounced  et  the  lowest 
heights  and  causes  substantial  ray  bending  when  9  is  close  to  O* 
(tangential  path) j  thus,  many  fine  steps  ere  required  initially. 

Modeling  is  Introduced  for  the  purpose  of  predicting  the  mean 
of  EXT  propagation  effects  and  the  limits  of  their  variability 
from  readily  available  climatological  data  bases.  In  order  to 
apply  such  e  tool  in  an  optimise  manner,  it  is  important  to  con¬ 
sider  the  uncertainties  in  the  spectroscopic  date  end  the 
limitations  imposed  by  the  numerical  calculation  procedure. 

D.  Typical  Examples 

Ttie  examples  presented  in  Figs.  1  to  7  almost  speak  for 
themselves.  Figure  1  ia  e  presentation  of  the  complete  rota¬ 
tional  water  vapor  spectrum  for  e  typical  sea- level  condition. 

Zt  wee  calculated  by  using  the  AFGL  Data  Tape  (10)  rather  than 
the  EHF  model  given  in  sections  8  end  C,  end  serves  to  demon¬ 
strate  the  very  strong  (over  6  orders  of  magnitude)  absorption 
due  to  water  vapor  in  the  far- infrared  that  allows  on  the  low 
frequency  end  just  a  few  window  ranges  for  t  ran  mission. 


The  continuum  absorption  If”  in  ppm  is  of  e  magnitude  end 
functional  form  that  is  similar  to  the  empirical  Gaut-Aei finite  in 
expression  (11) 

•  .  S.S  x  lo'5  p  •  t1'1  (v/30!  (27) 

which  wee  introduced  by  waters  (31)  and  hae  proven  useful  to 
other  workers  (32,108).  The  fer-wing  contributions  to  refrae- 
tivlty  of  all  rotational  HjO  lines  yield  with  Eq.  (12)  in  Eq.  (13) 
the  term  41.6  at2  (3,4,5).  The  contribution  of  tha  local  lines, 
i  •  46  to  74  to  thl*  value  is  in  ppm 

r[WMv-0)J.  ...IS.,5  (2«) 


A  more  refined  mode)l*g  result  is  depicted  in  Fig.  2.  Zt 
displays  at  higher  altitudes  (h  ■  16  km)  spectral  signatures  of 
the  tree*  geees  Oj,  CO,  N2O,  which  are  not  included  in  the 
described  model.  At  tropospheric  heights  (h  •  0  and  4  kxi) ,  only 
lines  of  HjO  end  02  (see  markings  at  bottom  of  figure)  are 
important.  The  figure  gives,  in  essence,  a  picture  of  atmo¬ 
spheric  molecular  absorption  in  the  millimeter  and  submillimeter 
vavelenr-.n  range.  Thee*  computations  by  Burch  end  Clough  (24) 
agree  well  with  reported  date  and  the  presented  model  (Eqs.  (27) 
and  (4))  of  the  water  vapor  continuum  spectrum  labeled  o^.  The 
attenuation  rate  spans  saven  orders  of  magnitude  over  the 
frequency  rang*  v  -  100  to  1000  GHz  and  altitudes  h  *  0  to  16  km 
where >99.5%  of  the  atmospheric  water  vapor  is  contained.  The 
transparency  in  the  window  ranges  W2  to  W5.  which  are  valleys 
between  the  absorption  line  peaks,  is  dominated  by  the  water 
vapor  continuum  a,.  Figure  2  affirms  the  dominant  role  that 


TABLE  2. 


Data  Base  for  H^O 


Spectral  Lines  in  Air  up  to  1000  GHz 


ID  (lowr 
quint,  no.) 
(10) 


Frequency  .  GHz 


ricunw  5.  Zenith  ett mustion  A  throw?*  the  firm t  30  km  of  thm 
U.S.  Standard  Atmosphere  197 4  for  dry  (PH  -  St).  moderate  (3 0%) , 
•nd  humid  (100%)  sir  over  thm  frequency  range,  v  •  JO  GHt  to  300 
GHt  (Syr).  Thm  rmla tivm  humidity  pH  w sm  attummd  to  bm  constant 
for  h  -  0  km  to  t  km  (19) .  Thm  mwpmrimmntml  dstm  srm  from  Pmf. 

40  snd  Tsblm  5. 


Frequency ,  GHz 

riGUPB  6.  Zmnith  phatm  distortion  A4  -  0  -  60  (tq.  (It) )  for 
thm  «<nc  conditions  mpmcifimd  in  Fig.  5.  Thm  frmqymne y- 
indmpmndmnt  doing  timm  duo  to  #0lff 0(h)1  in  (Kg.  (9)3, 

PH.%  5.000  50.000  100.00* 

H,ns  7.600  7.977  9.26* 


water  vapor  absorption  plays  in  atmospheric  millimeter  wavs 
propagation,  avan  at  nod a at  humidities  (Ml  -  46%  for  h  •  0) . 

A  norm  detailed  picture  of  Molecular  attenuation  appears  In 
tig-  J.  The  oxygen  microwave  spectrum  dominates  In  the  ranges 
SO  to  70  GHz  and  11S  to  123  GH*.  At  levels  cloae  to  the  surface, 
the  60-CMx  lines  are  merged  into  an  unstructured  band  shape,  the 
manimm  intensity  of  which  la  pressure-proportional,  until  the 
lines  separate  (h  >  IS  km) .  Above  IS  km,  the  shielding  effect 
breaks  down  and  radio  channels  with  up  to  400  MX*  bandwidth  can 
be  accoemdated  between  the  lines.  Zf  there  is  a  line  close  to 
the  frequency  of  interest*  isolated  line  behavior  takes  over.  In 
the  frequency-agile  applications,  it  is  possible  to  tune  to  a 
more  or  less  constant  shielding  factor  (for  example,  0.S  ds/km) 
over  the  height  range  h  •  0  to  *  JO  km.  Above  h  >  30  km,  tenon 
splitting  has  to  be  considered  in  the  frequency  intervals 
Uq  i  S  MKs  until  the  O2  lines  vanish  above  h  >  100  fcm  (19) . 


frequency.  GHz 

rictmz  3.  Attenuation  rate  set  thm  altitudes ,  h  •  0,  5,  io, 
JO  km  over  thm  frequency  ranqm,  v  -  49  Offs  to  125  GHt  ditplay/nq 
bond  mod  linm  structure  of  thm  omygmn  mieromsvm  spectrum. 


riGVPt  4.  pure  water  vapor  tmfraetimity  nQ  ♦  Dfu;  fOT 
a  •  1.705  APe  over  thm  frequency  range,  v  *  3 00  Offs  to  J 000  Gffs 
from  measurements  end  calculation  (67; . 


tt  la  encouraging  to  the  author  to  tee  experimentally  verified 
results  on  milllmetsr  wave  water  vapor  refractivlty  axpreased  by 
No  ♦  0(v),  an  example  of  which  is  prssented  in  Pig.  4.  After  many 
careful  measurements  of  microwave  refractivlty  in  the  fifties  snd 
sixties  (2, 1,5, 47-50. $6) ,  this  marks  s  new  start  in  a  heretofore 
inaccessible  frequency  rsnge.  In  comparison  with  conventional 
absorption  spectroscopy,  refractive  dispersion  measurements  are 
superior  for  absolute  intensity  studies  (15,57,63).  An  important 
rssult  of  Pig.  4  Is  the  fact  that  the  much  stronger  llnea  beyond 


0  »  m  iso  m  no  n 

Frequency*.  GHz 


riGVPS  7.  Zmruth  brightness  temperature  (atmospheric  noise) 

Tg  domnwollinq  to  ground  from  e  dry  (w  *  0)  and  a  noise  (v  -  2  cm, 
rg.  (33))  ear  mats  (31) . 


A  NEMESIS  F0«  MILLIMETER  WAVE  PROPAGATION 


1000  CHs  determine,  by  means  of  their  rv»  intensity,  seat 

of  tha  microwave  rafr activity  mq  |f««  Eq .  (.'*)). 

E.  Cumulative  Behavior 

A  standard  example  for  an  inh-wogeneous  meditmi  is  tha  one- vs y 
sanith  response  through  tha  U.S.  Standard  Atmosphere  197ft  (19). 
figures  5  and  6  display  tha  cumulative  attenuation  A  and  cusula- 
tiva  phase  dispersion  A5,  for  which  48  height  levels  up  to  h  •  10 
km  are  simwcd.  Except  for  the  vicinity  <ve  -  10  tdis)  of  tha 
Oj-HS  lines,  these  curves  represent  the  one-way  sanith  path 
behavior.  Three  humidity  profiles  D(h>  were  used  to  model 
W  *  $i  50,  and  100%  in  each  of  2ft  height  layers  between  h  •  0  and 
•  km.  The  m  drops  rapidly  below  1%  above  h  -  8  km.  The 
frequency  range  55  to  85  CM*  is  opaque  (A  >  30  d8)  for  any  system 
attesting  to  look  through.  The  transmittance,  Eq.  (7)  is 
measurable  when  the  atmosphere  is  somewhat  transparent  (A  <  30  dll, 
and  it  can  be  determined  either  by  the  absorption  of  a  sienal 
coming  from  the  outside  (sun  emission,  satellite  beacon)  or  by  the 
thermal  emission  originating  predominately  in  8  km  thick  layers, 
frequencies  around  57  and  61  CHs  yield  the  maximem  and  minimum 
phase  dispersion  for  the  60-Oiz  band,  which  is  almost  independent 
of  the  water  vapor  content.  Actually,  the  foreign-gas  broadening 
of  the  Oj-MS  by  water  vapor  (see  Eq.  (22))  reduces  the  attenuation 
slightly  with  increasing  vapor  pressure  (see  fig.  St  A(10O%)  - 
A  (5%) ) .  Comparison  of  ths  total  attenuation  A  with  measurements 
(see  Table  5)  yields  reasonably  good  agreement  when  the  empirical 
•ssjvptions  formulated  by  Eq.  (27)  art  used  in  the  calculation. 

The  amount  of  water  vapor  absorption  for  a  given  atmospheric 
condition  can  be  eetlnatcd,  to  some  extent,  if  the  surface  vapor 
concentration  p^  la  known. 

The  EXT  thermal  emission  by  H^O,  02,  03  and  CO  was  calculated 
by  uetera  by  using  Eqs.  (10),  (7)  and  (4)  and  ths  absorption 
coefficients  given  in  Eq.  (31).  Calculations  are  shown  in  fig.  7 
both  for  no  water  vapor  and  for  a  total  preci  pi  table  water  vapor 
of  2  cm.  The  integral  is  evaluated  for  h  -  0  to  60  km.  Clear-sky 
emission  varies  primarily  with  the  snount  of  water  vapor.  The 
calculations  ere  for  observations  in  the  zenith  direction  from 
ground  level.  Cosmic  radiation  of  3  K,  which  is  incident  on  the 
atmosphere  from  the  top,  has  been  added  to  Eq.  (10). 

III.  PHYSICAL  PROPERTIES  OF  ATMOSPHERIC  WATER  VAPOR 

The  clear  atmosphere  appears  in  the  EKP  radio  path  modeling 
scheme  only  aa  a  p-t-e  (Section  XX,  B)  parameter  system.  The 
real  atmosphere  is  an  enormous  theater  of  diverse  but  related 
and  incessant  activities  (33).  Water,  endlessly  chancing  its 
phasas,  is  ths  main  actor  on  this  stage.  On  a  global  scale,  1 

each  year  about  1  m  of  the  ocean  depth  is  evaporated,  whereas  I 

the  average  water  content  of  the  atmosphere  represents  e  depth  1 

of  roughly  2  cm.  Water  vapor  ie  rapidly  exchanged  causing  it  to 
be  *  patchy,  capricious  meditmi  with  parcels  (scale  sixes,  10  to 
100  m),  blobs  (<  0.5  a)  and  strata  <»  1  m)  of  molsturs  that 
probably  account  for  a  fair  share  of  the  scatter  end  inconsis¬ 
tencies  in  millimeter-wave  propaoetion  data  taken  from  field 
observations.  The  total  integrated  water  vapor  end  liquid  water 
in  a  vertical  colimn  can  be  detected  with  high  time  resolution 
U  10  s)  by  ground-based  redlometry  in  the  micr«#ave  range  (40). 

Xn  relating  millimeter -wove  propagation  and  atmospheric  water 
vapor,  one  hae  to  be  generally  content  to  make  predictions  with 
no  more  than  statistical  certainty. 

The  fact  that  the  three  phases  of  water  contribute  to  the 
veatj»«r  cycle  at  the  prevailing  conditions  (p-t-e)  is  most 
fundamentally  a  consequence  of  the  molecular  structure  of  H20. 
water  molecules  tend  to  associate  throuoh  hydrogen  bonds-*  having 
al  —**  ono- tenth  the  strength  of  e  molecular  bond.  Very 
st  - .  tcally,  in  the  solid  phase  four  H -bonds  form  a  rigid 
lavw.ee>  in  the  liquid  phase,  on  the  average,  two  H- bonds 
amalgamate  chainlike  links >  and  in  the  gaseous  phase,  a  chance 
exists  to  associate  singly  H- bonded  molecules  to  dimers  (21,43. 

44)  or  aggregate  into  clusters  of  preferred  sixes  (that  is,  10  to 
50  molecules)  (22.28)  under  the  Influence  of  ion-activity  (45). 

Xf  further  notice  is  taken  that  atmospheric  air  is  never  free 
from  invisible  particles  having  a  variety  of  origins,  chemical 
compositions,  sites,  end  affections  for  water  vapor,  then  the 
radio  engineer  will  almost  despair  at  the  prospects  of  ever 
putting  order  in  the  atmospheric  pendemoniias.  The  value  of  the 
interdisciplinary  workshop  on  Atmospheric  Mater  Vapor  for  pro¬ 
viding  guidance  with  this  task  must  be  stressed. 


p<h)  *  po  **cp(-h/2  fcnj 

end  the  most  realistic  measure  for  predictions  is  the  total 
precipltable  water  vapor  in  as  f 1 0  g/r3/>x*  -  l  cm). 


?(s)  do  *  0.2  J 


The  quantity  w  is  measured,  for  example  for  e  sanith  path,  with 
a  microwave  radiometer  (39,40).  The  variations  of  the  dry  air 
parameter  p(h)  end  t(h)  are  described  by  standard  height  profiles 
end  ettentlon  is  celled  to  the  model  of  the  U.S.  Standard 
Atmosphere  (for  example,  03)).  The  number  density  for  the  atmo¬ 
spheric  gases  with  constant  vohr*e  mixing  ratio  (Oq,  W2,  C02 * 
nobis  gases*  or  referred  to  as  dry  air)  follows  directly  from  the 
p-t  combinations.  Over  the  heioht  range  h  •  0  to  16  km,  the 
changes  ere  for  p  from  101  to  10  kPa  (1  kPa  -  10  mb)  and  for  t 
from  0.9  (333  K)  to  1.5  (200  E) .  Ozone  (Oj) .  besides  water  vapor, 
hae  e  variable  mixing  ratio  which  is  modeled  separately.  Examples 
of  0^  millimeter  wave  spectre  are  shown  in  Pigs.  2  and  7  as 
reported  by  others. 

Fractional  fluctuations  due  to  water  vapor  pat****  and  tur¬ 
bulence  ere  typically  in  the  ranges 

6p  •  ♦  80% 

6t  -  t  1%  (34) 

jp  •  i  O.U 

Variability  of  humidity  for  a  mdlatltude  location  with  a  medium 
of  7.4  9/m-*  is  diurnally  -  1.0.  seasonally  t  6.3,  and  locally 
S  1.2  ,/■*. 

2.  Molecular  Quantitiaa.  Ths  molecular  world  of  sea  level 
air  is  very  empty.  The  molecular  radius  r  T  1.8  x  10~4  urn. 
average  spacing  d.  and  average  distance  t  between  collisions 
are  in  the  radios  r  i  d  :  X  •  1  :  19  *  320  (33) -  most  of  the 
time,  e  molecule  is  unperturbed  by  neighboring  molecules.  The 
nuetoer  of  dry  air  molecules  per  unit  volume  tm~3)  i« 


(e.g. ,  p  ■  101,  t  •  1  i  ■  2.44  x  10*5). 

Water  vapor  ie  an  imperfect  gas.  From  therxxxlyramc  measure¬ 
ments.  it  is  known  that  there  are  slightly  more  H?o  molecules 
per  unit  volume  than  predicted  by  the  ideal  g an  law.  The  correc¬ 
tion  is  made  by  introducing  the  second  virial  coefficient  8U>. 
The  molecular  nurt>er  density  of  H20  for  a  given  vapor  concentra¬ 
tion  (in  g/m3) 


follows  from  the  relationship  in  m 

*y  .  tta.wvp)  »•*»  ♦  »iuf 1  *  J.3«  »  io”p 

(e. g. ,  0  •  10  >  My  >  3. 35  x  1023) . 

very  few  values  for  »(t)  have  been  reported  (44,69) t 


a.  Absolute  Humidity 

The  amount  of  water  vapor  (that  is,  absolute  humidity  o) , 
present  in  the  atmosphere  depends  upon*  (a)  evaporation  from 
‘surfaces i  (b)  transport  by  notions  on  various  scalss,  mainly 
through  the  troposphere*  (c)  condensation- forming  clouds  and 
fog  that  causes  precipitation  fallout.  The  radio  propagation 
engineer  measures  the  amount  of  water  vapor  by  means  of  an 
average  concentration  p,  which  varies  at  sea  level  between  the 
extremes  of  0.1  (dry,  winter,  polar)  and  60  g/n3  (wet.  summer . 
tropical).  The  height  distribution  is  approximated  from  a  known 
ground  level  Po  by 


t . 

0.9 

1 

1.2 

T,  K  .... 

333 

300 

2S0 

W-*V  . 

-2.8 

-1.9 

-0.9 

deviations  from 

ideal  behavior 

c  •  lH/K°  l» 

■ 

o 

• 

(30) 

small,  even  at 

saturation 

(see  Eq. 

(431) * 

t . 

0.9 

1 

1.2 

p  *  0,.  <>/»5 

. 

130 

25.5 

0.82 

t.  10°  .  .  . 

12.3 

1.6 

0.03 

The  deviation  c  led  to  the  postulation  of  a  dimer  molecule 
t*20)2.  The  molecular  structure  (43)  end  the  millimeter  wave 
spectrum  of  the  dimer  (6,8,69)  are  well  established*  however, 
its  number  density  under  tropospheric  conditions  is  not  known. 

Ths  dimer  number  density  ie  expected  to  depend  strongly  on 
temperature  since  the  hydrogen  bond  strength  is  rather  weak. 

Data  which  era  suspected  to  be  caused  by  e  dimeric  effect  can  be 
analyzed  for  their  percent  chenge-per-degree  or  for  their  power 
law  (tF>  dependence  to  allow  e  positive  identification* 


Jtho  electron-rich  end  of  O  in  the  polar  molecule  AjO  attracts 
an  electron-poor  9  end  of  a  neighbor  HjO. 


f.t.' 
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eV/molecule  kcal/mole 


Data 

%  Change/K 

>1.56 

-2.54 

-3.7« 

-6.22 

-12.10 


tond  strength  la  expressed  in  units  of  ev/mDlecuie  or  kcal/mole 
and  the  v«iu«  «*p*ct*d  for  a  hydrooen  bond  lias  between  3  and  $ 
kcal/mole.  Th*  dinar  number  daniity  in  was  proposed  to 
follow  an  Hprmiofl  (6) 


An  approx i nata  aquation  in  ■»•*  for  k(t)  waa  derived  f row  axprea- 
aiona  given  by  Boh lander  (69) 

k  :  l.l  .  to'J1  tV*  (401 

which  whan  combined  with  £4..  (3?)  .nd  (391  .  yi.ld.  in  .  * 

h„  =  2.0  «  ID1*  0 1  t5'6  <413 

D 

and  for  the  fractional  dinar  concentration 

HIJ/Hv  »  6.1  *  10“5  0  t5'6  M21 

(for  exanple,  0  •  20.  t  •  1  :  * _/H  *  l*2  x  10  . 


Other  estimates  on  tha  fractional  dinar  concent  rat  ion  Mt>/*v 
are  report  ad  in  Ref.  (21):  0.4  x  10" 3  for  0  •  20  assuming 
$.4  kcal/nolaj  and  in  Raf.  (44):  1.9  x  10"J  at  t  •  1  and  0.9S  x 
10“ 3  at  t  a  1.11,  when  extrapolated  fro*  saturated  water  vapor 
data  taken  over  tha  range  t  •  0.77  to  0.84.  Thaaa  ratios 
Mp/Hy  are  cloaa  to  values  of  C  defined  in  Eq.  OB). 

B.  Relative  Humidity 

Tha  eiaaplaa  of  a  taw  spheric  nil  line  ter  wave  propagation 
(Pigs.  1  to  7)  showed  clearly  tha  doeunant  influence  of  water 
vapor.  Sane  difficulties  in  predicting  this  Influence  obviously 
are  related  to  tha  fact  that  water  vapor  is  s  vapor  and  not  a 
gas.  The  first  consequence  of  tha  vapor  state  is  that  s  nsxisnss 
concentration  0,  cannot  be  exceeded.  A  balance  exists  between 
two  states  called  tha  saturation  point.  At  saturation,  the  rote 
at  which  aolacules  evaporate  fro«  a  plane  surface  (liquid  or 
solid)  equals  the  rate  of  lnconing  condensing  Molecules.  The 
saturation  concentration  over  water  for  at*o*rhcric  conditions  is 
fitted  by 

P  <  •  17.39  t6  10110-’-'"'’  (43) 

in  .nd  la  In  a  vary  rough  approximation  airily  in  ,/«7 

P,  :  J5  4  11  (44) 


Equation  (43)  is  progressed  into  tha  radio  path  Model  (Section  II), 
which  uses  RM(h)  information  to  calculate  with  t(h)  the  values  s 
or  o. 

Tha  saturation  concentre tion  defines  rslative  humid toty 

f  -  f/f%  -  0.04ptl?  «  0.29  «tl*  <_  1  (45) 

or,  in  tha  More  familiar  percent  notation,  RH  ■  lOOf  £  100%.  The 
wsxlMuM  vapor  concentration  varies  over  s  wide  range  ss  s  function 
of  temperature: 


Smi thorn an  Hydro metric  Table  (sea  (33)). 

^Triple  Point 

Tha  interdependences  between  RH,  w,  and  T  are  shown  in 
Pig.  8.  .  The  humidity  condition  RH  -  100%  is  s  delicate  balance 
point  for  phase  changes.  A  drop  in  temperature  by  1  X  decreases 
the  water  vapor  density  of  saturated  air  by  about  6%.  Energy  is 
released  during  condensation  whieh  was  stared  in  the  random 
motion  of  HjO  molecules.  One  cubic  meter  of  a  forming  cloud 
that  converts  1  g  of  vapor  into  the  liquid  phase  releasee  2.5  kJ. 

The  latent  heat  release  can  generate  electric  fields. 
Depending  upon  the  suddenness,  the  amount  of  cooling,  and  the 
water  vapor  supply,  mere  or  less  violent  updrafts  (such  as 
cumulonimbus  clouds  in  thunder storms )  feed  s  turbulent  air 
motion.  Small  aerosols  ara  loftsd  into  the  upper  region  of  e 


Water  Vapor  Concentration  p,  g/m^ 


tt^r- 


^  Precipirabie  Water  vapor  w,  cm 

riGVKt  8.  Tha  amount  of  precipitsble  water  vapor  *  in  a 
radio  path  of  length,  L  •  1  km  and  the  corresponding  homogeneous 
water  vapor  concentre tion  o,  both  as  a  function  of  relative 
humidity  Ml  for  varioua  temperature*  T.  The  temperature  depen¬ 
dence  of  saturation  mater  vapor  pressure  eM  and  maximum  vapor 
concentration  Pj  (Wd  m  100%)  (33)  is: 


r,  c 

e9,kPa 

p,.g/* 

-30 

0.052 

0.453 

-20 

0.125 

1.070 

-10 

0.286 

2.360 

0 

0.611 

4.840 

10 

1 . 227 

9.390 

20 

2.337 

17.270 

30 

4.243 

30.310 

40 

7.378 

51.020 

The  broken  lines  indicate  schematically  the  range  where  water 
uptake  by  aerosol  takes  place  (Fig.  9) . 

forming  cloud  while  the  larger,  heavier  ones  remain  suspended  st 
lower  levels.  The  snail  aerosols  carry  a  positive  charge,  the 
lower  levels  take  on  a  negative  charge.  The  charge  eeperatlon 
generates  high  electric  field  strengths  end  lightning  discharges 
occur  when  a  value  of  about  30  HV/n  is  reached.  Even  before, 
additional  ions  must  be  produced,  which  are  suspected  to  have 
catalytic  influence  on  the  formation  of  honowo locular  cluster 
(H2°)n  with  B  >  10  to  50  (28). 

C.  Submicron  Hydrometeors 

Atmospheric  air  is  never  free  from  invisible  particles 
(aerosols)  having  a  wide  variety  of  origins,  sixes,  and  chemical 
composition,  and  most  importantly,  having  the  ability  to  convert 
water  vapor  into  submicron  hydrometeors.  TO  account  for  aerosol 
activity,  one  must  follow  the  evolution  of  particle  sixe  spectra 


•** 

5 

if- 

x  r* 


>0.)  i  r 

ii  T 

j 

♦CN.  \ 

.  >V  y 
S,,.9»  \ 

O  *  •  \  Ciwer  v 


MOIST  Am 
a  *  0  am 


*.  i  m 


uj-  ■  **  i  «oo*x. 

I  nw-Jwium  in 


x\ 

'O’4  V*.\ 


T.  X 

t 

t-*7 

fit.  tq.  43 

?,•  M-  «4  | 

333.30 

0.900 

6.000 

130.350 

130.400 

350  Z 

300.00 

1.000 

1.000 

25.500 

25.500 

25  ; 

273. 16^ 

1.098 

0.203 

4.850 

4.840 

5  1 

250.00 

1.200 

0.045 

0.822 

0.822 

1 

m - «*9*OOto  *02% 

I  Smew,  irvmiibe  mm.  »  5*1* 


Lxx«J  mtne»  Commit 


\a%Y 


Slow 

SlrOM-CUMAM 

0.8 


I  40 

-  X  Momcuvt 


I  -ue»  j 
7-Moearaw  i 

3-Ct0u0bv»w  "i 


partidt  Rodius  log  Mnt 

FIGURE  9.  Schematic  site  distribution  of  number  density  ff 
for  equivalent  spherical  atmospheric  particles  from  molecular  to 
rain  drop  sixes,  r  -  1(T 4  to  10*  u*  (32.33.41). 
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and  th^  attendant  microphysical  processes  of  rasa  transfer  over 
many  orders  of  magnitude  la  sls«  (r  *  2  x  10“®  to  1  cm)  and 
number  density  <D  *  1053  to  1  m"3) .  This  problem  of  acala  It 
sketched  in  fig.  9  to  illustrate  the  extent  of  the  difficulties. 
The  amount  of  relevant  literature  dealing  with  aerosols  is  over¬ 
whelming  and  the  referenced  material  (33-451  it  probably  an 
incomplete  selection  in  the  search  for  nonmoleeular  nil  line  ter 
wave  absorbers  in  moist  air. 

Aerosol  models  start  from  dry  particles,  distributed  over 
three  distinct  site  ranges  (fig.  91  (35-3®) «  B,  nueleation  moder 
A,  accumulation  mode*  C.  coarse  particle  nodei  and  X,  hypotheti¬ 
cal  submicron  distribution  required  to  explain  EMA  (sea  Section 
XV) .  Each  mode  le  described  by  representative  value#  of  a 

log-normal  distribution  (23.43) .  A  major  process  for  the  forma¬ 
tion  of  the  n  mode  la  the  trace  aas-to-partlcle  conversion. 
Sulphur  dioxide,  for  example,  nucleates  with  water  vapor  to  form 
sulfuric  acid  primary  particles  in  Urge  msnbcrs  (as  high  aa 
1015  a~3) (  which  are  rapidly  (within  ms)  coagulated  by  Brownian 
notion  into  the  D  mode  (34).  This  mode  may  contaie  sizable  mbs 
concentrations,  0*(D|  >  20  ug/n3  { 3*0  ,  but  more  to  the  point,  it 
presente  an  enoreous  surface  area  S*  to  the  vapor  phase  and 
enables  the  aeroeol  to  respond  quickly  to  changes  in  relative 
humidity  RN.  Particle  number  density  m*  and  total  surface  area 
Sjt  for  1  me  of  water  distributed  in  spherical  droplets  of  radius 
r  withe 1  m3  are  related  as  follows 


r,  um 

D^»  m>iec/m3 

v 

N 

0.01 

1.4  >  1014 

3000 

A 

0.10 

3.4  1  1C11 

300 

C- 

1.00 

2.4  x  10® 

30 

The  mixture  of  gee  and  suspended  particles  is  celled  an  aerosol. 

Solution  droplets,  such  as  H2S04,  are  highly  hygroscopic) 
that  is,  they  adapt  their  size  by  taking  up  water  if  the  HD  in 
the  ad) lent  air  Is  larger  than  the  equilibrium  m  over  the  drop¬ 
let's  surface.  Xn  e  reverse  situation,  their  slse  will  shrink 
because  of  evaporation.  Other  hygroscopic  agents  are  salt 
particles  found  in  maritime  (NaCl)  as  well  as  in  urban  (for 
«x*«g>le,  (1*114)2904]  environments.  These  crystals  undergo  a 
sudden  phase  transition  to  become  solution  droplets  at  critical 
m  values  (m  •  76%  and  80%  for  NaCl  and  (WlgTySDg,  respectively). 
The  aqueous  particles,  or  so-celled  hydrometeors,  collect  in 
stable,  HD-dependent  sizes  predominantly  In  the  A -mode.  A  third 
type  of  aerosol  particle  is  relatively  passive  dust  grains  in 
sizss  cossisnsurate  with  the  A*  and  C  modes.  The  solid  matter 
accumulates  s  film  of  water  on  its  surface.  All  the  particle 
population  is  able  to  aggregate  water  (CD)  but  the  smell, 
hygroscopic  parts  of  the  CD  (CCD)  play  the  dominant  role. 

The  total  aeroeol  ness  concentration  in  tir  with  e  humidity 
of  0  -  1  g/m3  can  range  from  oh  -  1  (clean  air)  to  400  ug/m3 
(polluted).  At  m  •  50%,  an  average  mesa  loading  of  about  0.1 
mg/m3  has  been  deduced  from  data  collected  at  291  locations  in 
the  o.S.  (3*).  Values  between  0.5  (clean  air)  end  3  mg/m3 
(polluted)  ere  reported  for  C*  in  cloud  free  aerosol  (42). 

Water  vapor- to- liquid  conversion  becomes  effective  when  M 
exceeds  80%.  An  aerosol  population  can  triple,  even  quadruple 
its  dry  state  size,  end  thus  lead  to  more  then  e  hundred-fold 
increase  in  mass  concentration.  Since  the  average  size  stays 
below  1  um,  a  sere  ox  less  Invisible  cloud  exist*,  optically 
categorised  in  the  atmosphere  as  haze.  Optically  active  sizss 
0.3  um)  in  sufficient  number  ere  only  produced  If  the  limit 
RH  -  100%  is  exceeded,  in  that  case,  the  available  amount  of 
water  (for  example,  that  due  to  n  sudden  cooling  of  saturated  sir) 
n  shared  either  by  many  (>  10*)  active  particles  (CCD)  and  small 
tr  *  B  isa) ,  uniform  droplets  fora  as  fog  or  cloud,  or  by  very  few 
t«  10®)  and  large  (r  >  100  um),  destabllslng  droplets  that  pro¬ 
duce  rein.* 

The  water  uptake  of  aerosol  has  a  theoretical  foundation  in 
tha  relationship  between  the  ambient  RH  and  the  equilibrium 
radlur  r>*  a  particle.  Hanoi  provided  empirical  data  on  the 
hydro^ctrlc  activity  of  aerosol  (35).  He  defined  e  mess  growth 
factor  q(f  -  nt)  with  respect  to  the  dry  mass  concentration  Cj 
and  measured  in  the  laboratory  the  equilibrium  growth  behavior 
of  typical  air  samples.  Two  examples  are  depicted  in  Fig.  10. 

The  accretion  of  water  to  pert  idea  ia  substantial  when  Ml 
exceeds  80%  and  gives  them  almost  pure  water  properties)  that  is, 
p*(f  >  0.8)  -  a**  TVo  approximate  expressions  for  g(f)  can  be 
formulated  (35) 


if  differences  end  hystareees  in  the  growth  behavior  due  to  the 
chemical  makeup  ere  disregarded.  A  detailed  discussion  of  models 
for  growth  factors  earn  be  found  in  kefs.  38  end  3T. 

By  referring  to  radio  path  modeling,  it  is  possible,  at  this 
point,  to  formulate  e  path- in teerated  liquid  water  content  in  " 

w  •  j  ds  •  j  PJ»  gff)  ds  «  j  lpj/«l-f))  ds  (4f> 


AfhQSX  «MTF*  UST4SC 
20-  C 


WSm  ******  1  mt,  x 

ricu* t  10.  water  uptake  (mas#  growth  /actor  gff  -  xx;j  of 
two  representative  aerosol  samples  (maritime*  NaClt  Continental : 
,nh4> 2*°4  *  SlOj )  as  a  /unction  of  increasing  and  docreamtnq 
relative  humidity  XX  (35) . 


Tha  water  content  w  increases  at  the  expense  of  the  precipltable 
water  vapor  u  defined  by  Eq.  (33);  hence,  in  cm 

w(f)  -  w(f-0)  -  W  (50) 

The  amount  of  liquid  water  in  moist  air  (RM  <  100%)  is  difficult 
to  measure.  Several  indirect  methods  have  been  developed.  Three 
more  promising  techniques  ere 

1.  The  mess  concentration  r£  of  a  dried  air  semrle  1* 
determined  by  micro-weighing  end  available  (35)  empirical  growth 
factors  gif)  ara  applied.  For  example,  it  TO  •  99%,  the  values 
for  g  range  between  65  and  120  and  yield  N  -  0.013  to  0.024 
when  o£  -  0,2  mg/m3  over  a  path  length  L  •  1  tat. 

2.  A  passive,  radiometrie  microwave  method  measures  at  two 
(or  more)  dtfflE*u*t  frequencies  the  sky  noise  T*  (see  Eq.  (10) 
and  ftafs.  39  and  40).  The  briohtnesa  temperature  T#  is  coneerted 
into  zenith  attenuation  A  (dB)  by  means  of  the  spectroscopic  data 
base.  Gulrsud  at  el.  (40)  perfected  the  technique,  which  uses  an 
instrument  operating  at  20.6  and  31.65  GHz.  The  retrieval  algo¬ 
rithms  for  the  simultaneous  determination  of  water  vapor  and 
water  liquid  content  ere  adjusted  to  take  into  account  e  priori 
statistics  (that  is,  a  representative  radiosonde  profile).  Tor 
example,  the  climatology  of  Oklahoma  City  yielded  (40)  in  cm 

w  •  6.22  A (20)  -  2.72  A(31)  -  0.01 


V  •  1.30  A(31)  -  0.53  A (20)  -  0.001 

The  frequency  20.6  CDs  is  elose  to  the  22  GRz  vapor  line  (Table 
2)  and  responds  mainly  to  vapori  the  frequency  31.65  CDs  is 
remote  from  the  line  end  responds  to  the  dielectric  loss  of  water. 
This  lnatnisnt  provides  around-the-clock  monitoring  of  w  and  M. 

A  typical  record  over  7  consecutive  days  (Oklahoma  City,  April, 
h  •  360  m,  locking  at  smith)  displayed  for  each  24-tour  period 
the  following  extremes  (40) 1 


2.5  to  1 
1  to  2.6 
2.4  to  0.8 
1  to  2 
0.5  to  1.2 
0.8  to  1.3 
2.4  to  4.2 


w,«w 

0.01  to  0.02 

0.02  to  0.5 
0.02  to  2 
0.03  to  1.5 
0.1 
0.1 

0.1  to  1.5 


The  point  of  this  example  Is  that  tha  low  limits  of  «  usually 
excaad  tha  valuas  that  are  obtained  when  using  an  averaga  dry 
aerosol  mass  concentration  end  assuaed  g(f)-  factors.  The  high 
limits  of  w  are  most  likely  due  to  visible  clouds  since  their 
size  range  lr  •  7  to  100  urn)  contains  the  bulk  of  the  water. 
Unfort"*»%tely,  the  data  wars  not  accompanied  by  notes  providing 
information  or  meteorological  or  optical  conditions.  A  similar, 
satellite- bom#  radiometer  yields,  over  ocean  surfaces  with  well- 
defined  emission  temperatures,  the  same  information  on  w  and  N  on 
a  global  scale  (39). 

3.  The  measurement  of  optical  transmission  gives  insight 
into  the  average  aerosol  state.  Both  the  optical  (X  -  o.SS  um) 
attenuation  rate  and  the  visibility  v  •  20  dB/tx.  are  related 
to  the  liquid  veter  content  0*.  In  general.  Die  scattering 
theory  has  to  be  applied  „n  order  to  predict  values  of  a,  based 
on  available  size  distributions  and  the  complex  refractive  index 
of  the  particle  materiel  (23,33,35). 
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Tor  the  purpose  of  this  r«P«.  **  !•  eutfleient  to  byp**»  the 
elaborate,  lengthy  calculations  and  discus*  the  rain  features 
with  **t  help  Of  a  simple  empirical  expression  (21,37) 

•u  5  300  Pw  6  °r  V  *  °-°*7/Pw'6  <52) 

where  la  In  dl/km  and  V  is  In  km.  The  relationships  between 
Oy,  V,  and  the  path-averaged  water  content  5*  are  approximately t 


K 

.  g/m3,.  < 

0.001 

0.01 

0.05 

0.1 

0.5 

%  • 

,  dB/km.  < 

5 

19 

50 

75 

190 

V  , 

km  ...  > 

4 

i.ru.si 

0.4  (0.5) 

0.27 

(0.21  0.10  (0.0«| 

Clear 

Hass  .  RH 

| 

vl 

fog. 

m  >  100% 

•l  1  •  experimental  value 


ITT 

1.  A  molecular  approach  searching  for  water  polymer*  (MjO)n 
and  their  ".prctTS  in  the  atmosphere.  Sixes  of  n*2  idrmer)4 
16,8,21.69,31,86,09)  and  n*)  (cluster)  (21.26)  are  po*  sib  Hitler. 

2.  Liquid  water  uptake  by  sufmeron  aerosol  particles  under 
conditions  of  high  relative  humidity  (RH  *  B5»l  (9. 20.23.25. 30, 

35). 

Each  conjecture  la  supported  by  *on»  as  eel 2  as  contradicted  by 
ocher  bits  and  piacas  of  experimental  evidence. 

A.  Definition  of  EWA 

window  attenuation,  both  model  and  experimental,  ia  fitted 
to  expression  In  dB/km  of  the  fore 

a  •  0.182VW"  •  C  p*  ty  (v/30> *p  (56) 

where  a*  y»  a  are  the  proper  exponents  of  a  particular  absorption 
modal. 


Frequently,  base  and  fog  conditions  are  described  merely  by 
stating  the  measured  value  of  the  visibility  v*  which,  however, 
depends  on  the  definition  of  the  threshold  value  (20  dS  for  L  • 

1  km  in  Eq.  (52)).  Nilsson  perforated  detailed  Hie  calculations 
on  five  different  aerosol  models  with  the  result  that  visibility 
V  (16  db  threshold)  and  relative  humidity  m  are  related  directly 
(23) .  The  range  of  KH  for  these  models  is 


v  .  km 

so 

15 

5 

2 

RB. 

20  to  73 

84  to  95 

96  to  99 

98.7  to  99.7 

In  slavery,  the  salient  points  of  this  section  arat 

1.  The  physical  foundations  for  path-averaged  water  vapor  w 
(Eq.  33)  and  liquid  water  content  *  (Eq.  49)  were  traced  is 
measurable  quantities.  Both  are  expected  to  be  the  naia  source 
for  millimeter  wave  absorption  in  moist  air. 

2.  The  number  density  h  of  the  major  molecular  absorbers  was 
olven  for  Oj  in  Bq.  (35).  for  H2O  in  Eq.  (37).  and  for  (lljOj 

in  Bq.  (41). 

3.  The  wa tar-vapor- to-wat«r-droplet  conversion  in  atmospheric 
aerosol  and  the  associated  scale  problems  ware  discussed  to  aid 

in  the  coeprehansion  of  the  stages  that  a  water  molecule  has  to 
undergo  before  precipitating  out  of  the  atmosphere.  The  key  role 
of  relative  humidity  RH  (Eqs.  43  to  48)  in  this  process  was 
lxiderllned. 

4.  The  interdependence  between  the  optical  properties  (Eq. 

52)  of  moist  air  and  the  liquid  water  content  la  not  unique)  it 
depends  strongly  on  the  aerosol  sice  distribution. 


The  ATGL  compilation  lists  38  350  HjO  Unas  from  20  Ofx  to 
331  TU  (10)  ordered  in  seven  bands  (14) ,  of  which  1809  rota¬ 
tional  lines  up  to  w  «  13  THs  contribute  to  the  millimeter  wave 
continuum  (17).  A  fit  to  these  results  yielded  Eq.  (26),  which 
can  bo  reformulated  into  Bq.  (56).  Magnitude  end  exponents  for 
the  lino  continue  are  given  in  dB/km/g/n^/kFe  by 

C¥  •  ».«  «  10~* 

ud  (57) 

».  *  u 

Vs  » 

*»  *  1 

These  exponents  are  based  on  the  molecular  rvj-ber  density  (Eq. 

37)  displaying  the  exponents  x  •  1  and  y  •  1.  Vlth  Eqs.  (56) 
and  (S7),  it  follows  for  the  s«"ple  conditions  above  (v  «  300, 
p  •  100.  0  -  17.  t  -  1.023)  that  ly  -  2.75  dB/km. 

The  ditmr  spectrum  is  discussed  in  detail  in  Raft.  69  and 
109.  Absorption  should  follow  dependencies  given  by  Eq.  (57), 
but  modified  for  the  dimer  number  density  (Eq.  41).  The  response 
predicted  in  this  fashion  in  dB/km/ (g/*3) 2  is 

C_  -  (value  of  Table  4) 

0 

and  <58> 

■d-1 

*D** 


IV.  EXCESS  WATER  VAPOR  ABSORPTION— EWA 

Wat ar  vapor  absorption  (Eq.  14)  in  the  millimeter  weve 
windows  atams  from  the  fact  that  the  intensity  of  ths  local 
minima  ter  wave  lines  N{  ■  EST*  (Table  2)  is  insignificant  and 
that  two  other  absorption  terms,  and  BJJ.  t scorns  dominant. 

The  far -wing  response  n£  of  the  rotational  spectnv  of  H2O  beyond 
1000  CNs-  (Pig.  1)  is  estimated  by  applying  the  approximation 
Y<<v<<v0to  ths  lime  shape  r"  (Bq.  16)  and  rssults  in 

r-  •  nv/v’  <sd 

where  P*  is  in  units  of  1/GBs.  The  H3O  far-wing  continum  of  one 
line  in  ppm  follows  with  Eqs.  (24),  (25),  (S3),  and  (11)  as 


*_  -  values  of  Tabls  4 
0 

Aerosol  liquid  water  attenuation  la  obtained  from  fxiblished 
dielectric  data  on  bulk  water  (Pig.  13)  yielding  in  dB/km  in  the 
Rayleigh  approximation  (c*  •  (n*)2  ♦  (ft*)2,  C*  *  2nin£)  (25.30, 
32). 

oA  *  0.»2vp(f  tVIK  »  2) “  »  CC“)  2|  (59) 

Another  method  of  calculating  the  aerosol  liquid  water  attenua¬ 
tion  assumes  that  the  medium  has  a  refractive  index  n£  *  1,  and 
that  the  attenuation  in  bulk  water  ia  increased  due  to  wavelength 
shortening  <lv  *  *©/*»*)  •  This  allows  to  formulate 

aA  *  (o ynM  (W/L)  (60) 


«;  ■  jt;  •  cta.so  ppt3*1  ♦  pa  taaw 


(54) 


where 

Ct  •  0.1,4  bj  bj/v^  (55) 

Per  moist  sir,  the  o2-term  due  to  self -broadening  ia  always 
smaller  (<  20%)  than  the  foreign  gaa -broadening  term  po.  The 
far-wing  contribution  of  the  strongest  M20  line  (Table  2i 
vQ  •  2774  GBs,  •  5.14  x  10*^®)  is  for  an  atmospheric  condition 
described  by  p  •  100  kPa,  t  •  1.023  (20*C).  0*  -  17.27  g/m3,  at 
the  frequency  v  •  300  GHs,  •  4.77  x  10-4  ppm  or  a,  ■  0.026 
dB/km,  to  which  the  c^contr ibution  is  11%.  The  line  canter 
attenuation  rata,  in  comparison,  ie  3*  ■  7.40  x  10*  dB/km.  Most 
certainly,  line  shape  theory  (Eq.  16)  is  overtaxed  when  it  ia 
applied  to  predict  relative  intensities  in  the  far-wing  over 
seven  orders  of  magnitude. 

An  unspecified  tarn  N£  was  added  to  Eq.  (14)  to  account  for 
discrepancies  between  predictions  based  on  summing  tty-terms  (54) 
and  measurements.  The  only  certainty  in  the  conflicting  evidence 
for  Nfl  is  its  correlation  with  atmospheric  humidity.  The  abbre¬ 
viation  ewa  (excess  water  vapor  absorption)  is  used  to  describe 
Ny.  For  simplicity's  sake,  different,  mostly  exponential  tempera¬ 
ture  functions  have  been  reduced  to  a  power  law  t?  (Eq.  31).  In 
this  section,  the  exponents  of  physical  (D*.  ty)  and  frequency 
(v*)  dependences  are  formulated  for  absorption  models  and  cree¬ 
ps  red  with  those  for  experimental  data.  This  procedure  appears 
to  be  one  practical  way  of  identifying  a  specific  absorption 
mechanic,  especially  since  EWA  observations  are  not  over- 
whelmingly  consistent. 


Both  swthods  have  been  applied  to  the  latest  dislectric  data  on 
water  (106)  producing  somewhat  different  values  (see  Table  4). 
Frequency  and  temperature  dependencies,  when  approximated  using 
Eq.  (60) ,  follow  In  di/smi  from  the  rough  data  fit 

o yn^  •  24  (V/J00)°'&s  t*‘ 

Mow,  If  the  liquid  water  content  Is  estimated  by  Eqa.  (49) ,  (47) , 
and  (45) ,  w  *  6  x  10“3  L  p2  t34  oj,  one  obtains  for  the  expected 
aerosol  response  in  the  form  of  Eq.  (56)  that 


(61) 


In  sunsary ,  millimeter  wave  window  attenuation  a  might  very 
wall  be  a  combination  of  up  to  five  different  contributions  1 


a  ■  (local  lines)  ♦  a^tfar -wings)  ♦  ♦  <»0  ♦  (cluster?) 

(62) 


A  parses trie  study  of  water  vapor  concentration  (p>  and  tempera¬ 
ture  (t)  dependencies  could,  in  principle,  reveal  the  following 
behavior  in  the  attenuation  ratei 


*jont  already  in  2919  had  considered  a  water  dimer  to  explain 
anomalous  dielectric  water  vapor  results  (72). 


-All- 


Two  schools  of  thought  have  evolved  to  explain  Mi 
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Multiplication  factor  with  respect 
to  300  K  data 

-j  “i  n? 


•c 

,/V 

(g/M  ) 

Eq.  (57) 

Eq.  (58) 

Eq.  (61) 
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Data  obeying  x  *  1  and  y  *  3  but  displaying  magnitudes  different 
fm  Cy  <Eq.  57)  can  b«  interpreted  as  failure  of  tha  line  shape 
function  (Eq.  16)  to  predict  far -wing  intensities.  Discrepancies 
of  this  nature  will  be  nost  pronounced  for  dispersion  intensities 
(Eq.  151 .  winq  data  of  0(v)  drop  wore  gradually  <«  l/\>)  with 
decreasing  frequency,  and  actually  blend  into  the  well-known 
refractivity  *«0  (Eq.  12)  for  v  <  100  CHs.  as  sssn  in  Piq.  4. 

Excess  Water  Vapor  (EUA):  Evidence  from  Laboratory  Data 

Laboratory  Measurements  play  an  important  role  in  verifying 
modeling  schemes  for  EHT  properties  of  moist  air.  Generally,  all 
tha  spectroscopic  parameters  (>  200)  entering  into  Eq.  (1)  should 
be  deduced  from  absolute  intensity  measurements  under  well 
controlled  v-p-t-o  conditions.  By  a  judicious  choice  of  the 
experimental  variables,  it  is  possible  to  investigate  nost 
parameters  separately. 

Numerous  millimeter  wave  and  infrared  studies  of  water  vapor 
and  aoiet  air  have  been  reported.  References  46  to  72  ars 
selected  for  their  bearing  on  the  EWA  problem.  Evidence  for  EVA 
from  these  efforts  is  summarized  in  Table  3.  Zt  is  not  limited 
to  absorption,  but  also  shows  up  in  dispersion  spectra  D(v)  (63) 
and  in  refractivity  No-studies  (41,49. 52. 56,63) .  The  experiments 
are  performed  by  various  technique*.  The  radio  path  is  simulated 
in  an  enclosure  either  for  a  single- transit  or  a  multiple- 
reflection  (resonator)  passage.  Detection  sensitivity  increases 
with  path  length.  The  transmitted  energy  can  be  a  single 
frequency,  a  frequency  pair  for  differential  measurements,  or  a 
broadband  (Fourier  transform)  signal.  Main  variables  for  an 
expariment  are  either  frequency  or  pressure,  the  latter  being 


preferable  for  EWA  studies.  Gas  mixture  control  and  vacuum 
reference  are  laboratory  advant  for  ah*  Mute  lnt*-n**ity 
studies. 

water  vapor  is  recognized  to  b»  i  m'dlum  that  in  difficult 
to  control,  even  in  the  laboratory,  due  to  its  attractive  force 
toward  surfaces.  The  surface  area  of  the  laboratory  enclosure 
.vplaces,  in  a  way,  the  micro-surface  of  an  atmospheric  aerosol 
population.  Water  molecules  do  not  ordinarily  aggregate 
spontaneously,  but  water  vapor  becomes  liquid  water  when  wet table 
surfaces  are  present  to  retain  the  impinging  molecules.  A 
threshold  of  >  2  x  10"*  cm"*  the  surface -to- volume  ratio  is 
sufficient  to  form  a  continuous  interface  between  vapor  and 
liquid  06),  This  value  is  always  exceeded  in  atmospheric  air 
(tee  Eq.  46)  as  well  as  in  a  laboratory  test  chamber.  An 
experimenter  can  select  surface  rater la is  that,  to  a  certain 
extent,  peaaivate  the  attraction  for  water  vapor.  A  systematic 
study  was  made  of  various  surface  coatings  applied  to  an  electro- 
polished  stainless  steel  (SS  304)  cavity  (3440  em*  and  1266  cm2, 
S/v  m  0.37  cm"l)  evacuated  for  >  24  hours  to  10"*  torr  and 
subjected  to  pure  water  vapor,  eQ  •  2.40  kPa  at  23.0*c.  Results 
on  tha  ralativa  amount  e/t0  of  vater-uptaka  by  the  walla  and  the 
time  response  of  the  absorption  ('■'*$)  end  desorption 
process  ars  as  follows: 


Coating 

*o.b 

\a 

T.f Ion  rono  (OuPmitl, 

% 

6 

& 

’V 

1  coat 

-1.00 

20 

180 

500 

SS  304,  elect ropoliahed 

-1.50 

30 

150 

400 

HMDS*  silanizinq 

-1.15 

27 

190 

750 

Parylene  C  (Union  Carbide) 

-1.60 

32 

140 

580 

Silicone  SR240  (GE) 

-2.10 

42 

20C 

550 

% 

Teflon  FEP120,  2  coats 

-2.90 

58 

230 

5500 

Clear  lacquer 

-4.00 

80 

300 

7200 

60  GHz  spectrometer  cell, 
untreated  (63) 

-19.  Sill 
(S/V  -  1. 

102 

.28  cm"1) 

3  hrs 

5  hrs 

*He»anethylldlailaxane  l ICHj) 

^Number  of  molecular  layers 

Tha  need  for  liniexi  S/V -ratios  and  judicious  matsrial  selection 
for  laboratory  enclosures  is  clearly  evident  upon  comparing  the 
performance  of  on*  of  the  typical  spectrometer  cells  with  the 
preceding  test.  Water  vapor  surface  effects  have  been  recognised 
(52) i  yet.  In  many  cases  (46,49,56),  they  were  excused  as  a 
possible  source  of  error  for  the  reported  data,  work  close  to 
saturation  requires  a  circulating  gas-handling  system  with 
controlled  mixing  and  continuously  monitored  ffl  levels  (107). 
Additional  sources  of  error  are  disturbances  in  thermodynamic 


TABLE  3 .  Summary  of  Laboratory  Studies  of  Continuum  H^0  Absorption 
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in 


equilibrium:  (a)  the  vapor  heata  up  vftan  injected  Into  an 
evacuated  cavity,  (b)  the  adsorption  process  releases  heat  and 
vice  versa,  (c)  the  aablent  gas  temperature  is  lowered  during 
purp-down.  All  laboratory  results  on  spectroscopic  data 
(W0,  D.  H")  of  water  vapor  should  be  seen  in  light  of  these 
cements.  Some  selected  examples  displaying  EVA  behavior  ere 
discussed  in  the  following. 

1.  Moist  Air  Studies.  Llewellyn* Jones  et  el.  (66) 

Invest igeted  the  temperature  dependence  a l  the  frequency 
v  -  213  ghs  by  studying  moist  nitrogen  (p  -  93  kPs)  over  the 
range  T  -  270  -  320  K.  Hie  attenuation  (Eg.  62)  for  these 
conditions  can  be  described  by  the  empirical  expression  in  dB/toe 

a  •  o.llpt 2  .  0.01fi2t4  *  0.03f2t26  («1) 

Close  to  saturation,  the  following  valuea  result  from  Eg.  (63) i 
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2.  Purs  *at*r  Vapor  studies.  In  the  caae  of  self -broadening 
an  exponent  x  ■  2  (Eg.  54)  is  expected.  Mrowinski  (56)  observed 

a,  >  )  x  10”4  P2  dl/km  at  V0  -  22.235  GHz  and  T  ■  39*C,  in 

addition  to  the  line  absorption  <S|.  Measurements  by  Lie  be  (63) 
et  30.6  CMs  and  61.2  GKx.  T  -  300  ana  325  K  yielded  in  dB/km 

a  *  4  *  10*3^t2(v/60)2  ♦  6  x  10-Vt10(v/«0>2*5  (64) 

‘Hie  same  experiment  gave  for  the  refractive  dispersion  (ppm) , 

0  •  N* (61.2)  -  N* (30.6) ,  the  result 

D  •  2  x  10’V8  ♦  1.2  x  10-Vt“  {65) 

wn_  .  is  an  average  of  the  type  of  result  exhibited  in  Pig.  11. 

The  dispersion  response  resembles  the  water  uptake  curves  (see 
Pig.  10)  published  by  Manel  (35)  and  provided  the  impetus  for 
the  dlscuasion  presented  in  Section  IZX.C. 

Boh  lander  experimented  extensively  in  the  100  GHs  to  1000  CHS 
range  end  deduced  from  attenuation  data,  including  results  from 
other  investigators,  the  component  o*/P2  shown  in  Pig.  12  and 
listed  in  Table  4  (69.  109).  He  also  calculated  (by  theoretical 


PXCURf  21.  Diapmraion  o  of  pure  water  vapor  measured  with  a 
spectrometer  cell  between  the  frequencies  61.2  and  10.6  GRa 
163).  A  strong  anomalous  component  is  displayed  exhibiting  a 
condensation/evaporation  hysteresis  typical  for  ms  ter  uptake 
activity  (see  Pig.  10). 


means)  the  envelope  of  e  rotational  dimer  spectrum,  a  fractional 
dimer  concentration  (Eq.  42)  of  10” 3  being  assisted.  The  data 
fit  between  EWA  end  the  dimer  attenuation  (e/x  in  Table  4)  is 
poor j  the  spectral  shapes  do  not  match.  The  fit  is  much  improved 
if  the  aerosol  attenuation  rate  (Egs.  59  and  49)  in  dB/km 

°A  * 

is  used  in  the  comparison  (s/x  in  Table  4) .  Equation  (66) 
assumes  pure  liquid  water  droplets  of  submicron  site.  Actually, 
based  on  this  fit,  it  is  possible  to  predict  a  liquid  water  con¬ 
centration  0*  that  is  needed  to  reproduce  the  CMA  data  of  Table 
4t  that  is.  in  g/m3 

•  1.1  x  10*5  pJ  (*■») 

At  this  point,  the  reader  might  recollect  the  arguments  that  have 
baen  brought  forward  to  formulate  Egs.  <4T) ,  (49).  (59)  - 
(61).  Millimeter  wave  attenuation  by  submicron  hydrometeors  is 
derived  from  Mia’s  scattering  equations  in  the  Rayleigh  approxi¬ 
mation  (33).  The  refractive  index  end  attenuation  rate  of  bulk 
water  are  presented  in  Pig.  13  (102-106).  The  bulk  water  attenua¬ 
tion  a*  la  higher  because  of  refractive  wavelength  shortening  in 
the  mediia*  hence,  the  division  by  n^  is  applied  for  the  gaseous 
aerosol  mediia  (a*  •  1). 


TABLE  4 .  The  EWA  Component  Ot^/p  Reported  in  Refs.  69  and  109  Compared  with  the  Water  Dimer 
Spectrum  Up/?2  Calculated  by  Bohlander  (69)  and  with  the  Aerosol  Liquid  Water  Attenuation  (a,b). 
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rictmt  12.  Suataary  of  reported  frequency  dependence  for 
excess  water  repot  ebaorption.  conpared  with  e  rotational  dimer 
band  a  pact  rum  (*&/*„  -  10~3)  (69)  and  the  absorption  a  pact  rim 
dj,  of  liquid  water  in  aubaucron  hydrometeora  (see  Table  4). 


J  30  300  GHz  3  30  THx  300 


Frequency w 

FlGVpt  13.  Complex  refractive  index  n*  Mnd  hulk  attenuation 
rate  ^  of  water  aa  a  function  of  frequency,  v  •  3  GHs  to 
300  Ttta  (composite  from  kef  a.  102  to  106). 

Tile  refractivity  N0  -  40. 4p  ♦  «N0,  when  Measured  at  40"C  And 
21.4  **u  and  23.6  GHz.  (49.56,61)  exhibited  -also  the  anomalies 
^o« 

RH.r  <  05  01.0  94.0  97.0  98.0  99.0 

6N0,  ppm  0  0.3  0.0  1.5  1.8  2.5 

C.  EWA  Evidence  from  Field  Measurements 

Millimeter  wave  field  measurements  are  carried  out  in  three 
ways  <37-101)« 

1.  attenuation  rata  a  <dB/km)  for  horizontal,  line-of -sight 
pa the' 

2.  total  zenith  attenuation  A(dB>»  and 

J.  sky  noise  Tg(K) ,  often  employed  to  infer  the  A  value. 

Ho»t  of  the  me  as or wants  are  performed  at  single  frequencies t 
a  few  were  carried  out  in  a  broadband  node  (Fourier  transform 
technique)  (86, 90*93, 98*101) .  Several  difficulties  plague  field 
observations! 

1.  absolute  calibration! 

2.  large  scatter  <±  lot  to  l  30t)  in  data  due  to  unspecified 
weather  along  the  path' 

3.  scarceness  of  data  at  high  Hiaiidltias  (RH  >  90t)  1 

4.  lack  of  simoltaneoua  recordings  of  the  integrated  water 
vapor  w  (tq.  33) > 

5.  absence  of  date  on  integrated  liquid  water  W  (Eq.  49) 
and  visibility  V  (Eq.  52)  in  cloud-free  airs  and 

6.  difficulty  of  fitting  data  empirically  to  surface-based 
meteorological  variables. 

Horizontal  and  zenith  path  data  are  suinarlzed  in  Table  5  and 
compared  with  model  calculations  by  using  the  u.s.  Standard 
Atmosphere  (19).  Aa  expected*  the  window  (Wl  to  HS)  absorption 
increases  with  the  water  vapor  concentration  0  or  path- Integra tad 
water  vapor  v.  it  was  observed  earlier  that  measurements  in  the 
100  GHz  to  117  GHz  range  revealed  considerably  higher  values  of 
water  vapor  absorption  than  were  predicted  by  R20  line  shape 
theory  (73). 

Generally,  the  absorption  A  is  divided  into  a  dry  contribu¬ 
tion  (Aj)  and  a  wet  term  (A*) .  The  dry  term  is  caused  by  oxygen 
absorption  originating  from  a  well-known  spectroscopic  date  base 
(Section  II) 1  hence,  la  *1  to  4).  it  can  serve  as  a  check  value. 


The  wet  term  determines  the  transr-arency  in  all  window  ranges. 
Transparency  Is  tightly  coupled  to  humidity  Ithat  is,  0,  w,  end 
W(RH)J  causing  various  propagation  limitation*  (for  example, 
usable  range  L  or  mlmcne*  elevation  angle  0)  with  in<~  resting 
amounts  of  water  vapor  in  the  rath  volume.  A  typical  case  esists 
when  weather  conditions  change  from  clear  and  dry  to  cloudy  and 
wet.  The  water  vapor  absorption  problem  fades  rapidly  away  above 
h  >  3  km.  It  has  become  a  standard  practice  to  separate  the  wet 
term  into  two  components!  one  proportional  with  0  (monomer),  the 
oU»er  with  a  0*  (dimer)  dependence.  The  first  fisld  results  to 
support  ths  dimer  hypothesis  were  taken  at  220  GHz  (76). 

Comparing  these  data  with  other  available  data  in  Table  5*  one 
notices  that  in  about  half  the  cases  a  squared  (z  -  2)  water 
vapor  dependence  improves  the  fit. 

A  horizontal  path  that  operated  at  182.9  CHz*  close  to  the 
183  GHz  H2O  lino  (Table  2)  *  exhibited  differences  in  the  water 
vapor  proportionality  depending  upon  clear  or  cloudy  sky  condi¬ 
tions-  The  difference,  3.90  compered  with  4.20  dB/kat,  wee 
reconciled  when  a  temperature  dependence  of  t10  was  assumed, 
which  is  about  four  tints  the  t  dependence  for  the  183  GHz  lino 
(Eq.  29*  t2*3)  (77).  The  same  experiment  produced  at  171  GHz  a 

relationship  between  attenuation  rata  P  and  concentration  that 
becomes  increasingly  nonlinear  when  0  >  10  g/w3.  Measurements  in 
MS  (330  CHs  to  360  CHz)  required  a  linear  (x  ■  1*  y  •  2.3)  and 
squared  (z  *»  2,  y  ■  11)  water  vapor  term  to  fit  the  data  (89). 
Condensation  effects  were  evident  in  recent  K4  end  W5  data  (101). 
For  the  first  time,  the  importance  of  relative  humidity  RH  was 
recognized.  An  EWA  contribution  with  z  *  2  end  y  *  16  to  30 
can  be  Isolated  and  an  increase  of  a*  with  beginning  fog  condi¬ 
tions  wms  measured.  Predicted  attenuation  rates  In  fog  are  in 
the  range  3  to  15  dB/km  for  v  -  100  to  300  GHz,  o*  *  1  g/e  and 
T  ■  4"C  (25).  An  independent  neasurerent  of  water  vapor  w  end 
liquid  W  contributions  would  be  desirable. 

Cumul  a  tire  data  taken  through  the  total  air  mantle  are  more 
difficult  to  interpret  for  their  phenomenological  origins.  In 
most  cases,  zenith  attenuation  A  is  correlated  with  the  surface 
water  vapor  concentration  CQ.  Simultaneous  measurements  of  the 
integrated  vapor  w,  either  by  microwave  (see  Eq.  51)  or  infrared 
techniques  (93)  arm  rare.  The  maximus  attenuation  detectable  in 
the  EKF  range  is  A  $  30  dB.  The  dry  term  Ag  of  a  slanted  radio 
path  follows  the  secant  law  Ad<90*)/sin9,  even  for  low  angles 
(9  J  10")  1  the  wet  term  A*  increases  very  rapidly  below  9  <  10* 
(91).  A  tangential  path  (9  -  0)  traverses  about  38  times  the  air 
mass  at  zenith  but  perhaps  100  times  the  water  vapor  content. 

It  seems  certain  that  A,  is  a  continua  spectrum  (79.93,98, 
108)  and  that  earlier  observations  of  spectral  dimer  features 
(81,86,90,64)  were  instrumental  effects.  To  quote  from  Ref.  93: 

* (EWA)  shows  no  spectral  features  snd  scales  with  w."  The 
frequency  dependence  of  the  empirical  continues*  spectrum 
(Eqs.  28  and  4)  was  s  •  2.  A  fit  of  the  Ay  term  between  15  GHs 
and  230  GHz  (Table  5)  yields  in  dB  approximately* 

*.■»  0.013  po(w/JOI1-7  (M) 

and  the  frequency  dependence  drops  to  s  •  1.22  for  window  data 
taken  between  140  GHs  and  300  GHs  (93). 

Field  data  mostly  serve  the  practical  purpose  of  establishing 
an  operational  data  base.  Agreement  with  ndel  calculations 
(Section  IX)  la  fair,  when  it  is  considered  that  the  model  stakes 
use  of  Eq.  (27)  to  describe  the  water  vapor  continuum.  Either 
integrated  vapor  w  or  surface  concentration  0  ere  useful  predic¬ 
tors  of  millimeter  wave  window  ebaorption  (Table  5) . 

D.  EWA  Discussion 

A  water  vapor  continue  spectrw  dominates  the  transmission 
behavior  of  atmospheric  window  ranges  (Wl  to  W5«  Figs.  I  to  3). 
Three  absorption  mechanisms  have  been  parameter ised  to  aid  in 
identifying  the  respective  relative  contributions  to  the  con¬ 
tinuum.  in  the  course  of  this  paper,  the  following  picture 
evolved: 
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The  empirical  expressions  ttqs.  62  to  68)  merely  organise 
experimental  datai  they  are  not  found  to  be  entirely  satisfactory 
to  vnlcualy  support  even  one  particular  absorption  scheme.  At 
present  it  is  difficult  to  relate  EUR  contributions  quantitatively 
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to  observed  variables,  Rvsults  vary  areally  fro*  on*  experiment 
to  the  neat,  even  under  controlled  laboratory  conditions.  EVA 
la  Moat  prevalent  when  AH  >  90%,  however,  the  best  type  of 
Measurement  haa  not  yet  been  undertaken. 

There  are  objections  to  both  interpretation*  of  EVA.  Unclear 
problems  with  respect  to  the  dimer  am  (a)  no  direct  physical 
evidence  in  the  atmosphere i  <bl  no  match  to  the  frequency- 
envelope  of  a  proposed  apeetrum  (Table  4),  and  (e)  no  unique 
•etch  to  the  teiaperature  dependence.  A  promising  EVA  hypothesis 
is,  et  present,  the  vapor-to-liquid  conversion  of  water  by 
aerosol  particles  under  conditions  vhei#  RM  ^  100%.  Her*  the 
objections  lie  in  the  fact  that  invisible  or  hare  cloud#  (see 
X -distribution  in  Fig.  9)  amst  exist  with  liquid  water  concentra¬ 
tions  close  to  those  neasured  in  fogs  or  clouds  (for  example, 

AM  *  96%,  T  -  10*C,  0  ■  9,  and  0*  •  0.1  g/*3,  by  using  Eq.  (67)1, 
Oelogne  raises  the  point  that  for  the  case  of  an  inhomogeneous 
distribution  of  liquid  aerosol  wster,  the  visibility  V  is 
Improved  due  to  the  0°' “-dependence  (Eq.  52)  whereas  millimeter 
waves  respond  just  to  the  total  water  v  (Eq.  49i  Ref.  10). 

which  explanation  for  EVA  has  the  most  merit  esn  only  be 
answered  by  additional,  completely  controlled  experiments  per¬ 
formed  in  the  laboratory. 

V.  CONCLUDING  REMARKS 

The  role  that  atmospheric  water  vapor  plays  in  millimeter 
wave  propagation  was  traced  by  means  of  data  from  modeling, 
laboratory,  and  field  studies.  The  aasisnptions  and  approxima¬ 
tions  made,  the  reasoning  used,  the  relative  importance  of 
various  parameters,  the  limitation*  of  available  data,  and  the 
separation  of  water  vapor  and  liquid  water  effects  have  been 
addressed,  mlecular  absorption  due  to  oxygen  is  prominent  in 
the  45  GHs  to  125  GHz  range  (Fig.  3) »  water  vapor  dominates  at 
higher  frequencies,  actually  up  to  30  TKx  (Figs.  1  and  2). 
Considering  only  molecular  absorbers  (H20  and  02)  for  a  cloudless 
atmosphere  leads  Invariably  to  discrepancies  between  predieted 
and  measured  attenuation  rates.  The  problem  is  most  apparent  in 
the  ERF  window  ranges.  A  horizontal  path  at  sea  level  can 
experience,  at  220  GHz  for  example,  the  following  varying 
attenuation  rates  (32) t 

Clear  air  1.6  to  11.2  dB/ka 

Fog  0.4  to  4.7  dBAm 

Rain  (<  10  am/Hr)  1  to  7  dB/km 

These  values  demonstrate  the  relative  isportance  of  understanding 
the  clear  air  problem  first.  Clear  air  attenuation  la  related 
to  atmospheric  water  vapor  content)  unfortunately,  more  than 
two-thirds  of  this  contribution  ia  described  only  by  empirical 
formulas  that  lack  both  physical  insight  and  general  appli¬ 
cability. 

•  A  reasonably  concise  model  for  the  molecular  absorption  was 

presented  in  Section  IX  and  then  applied  to  identify  the  magni¬ 
tude  and  frequency  and  physical  dependencies  of  the  non-HjO 
contribution  to  the  available  data  body  reported  for  laboratory 
and  field  experiswrnts.  The  contribution  is  called  EVA  (excess 
water  vapor  absorption)  and  It  was  found  the  EW»  exhibits  in  the 
millimeter  wave  range  a  continuum  spectrum  similar  in  frequency 
response  to  the  dielectric  lose  spectrum  of  wster.  The  search 
for  liquid  water  in  clear  air  led  to  the  water-uptaka  phenomenon 
of  aerosol  particles,  which  grow  rapidly  in  moist  air  with 
increasing  AM,  This  problem  haa  been  addressed  in  the  infrared 
(9,20,23)  and  deserves  further  study  with  respect  to  influences 
on  millimeter  waves. 

Moat  of  the  evidence  found  in  parameterizing  EVA  data 
supports  a  condensation  phenomenon  driven  by  relative  humidity, 
AM:  (a)  the  nonlinear  pressure  dependence  (aerosol  growth 
function),  lb)  the  strong  negative  temperature  dependence ,  (c) 
hysteresis  effects  when  cycling  the  AM,  (d)  the  governing 
variable  is  relative  (AH)  and  not  absolute  (P)  htaidlty,  and 
(a)  the  failure  of  molecular  approaches  hypothesizing  weakly 
bonded  dimers  but  neither  matching  tbs  spectral  response  nor  the 
temperature  dependence  and  mssber  concentration  necessary. 

From  s  practical  point  of  view,  one  ia  looking  for  invisible 
liquid  water  eon  tent  ia  the  range  V  «  0.02  to  0.2  am  when 
AM  i  100%,  which  are  normally  known  to  ealat  under  haze,  cloud 
and  fog  conditions.  But  clouds  and  fog  imply  only  larger, 
optically  active  (r  >  L  um)  particle  aizee.  Should  the  condensa¬ 
tion  hypothesis  have  merit,  then  CMA  can  serve  as  a  tracer  to 
the  otherwise  not  directly  accessible  world  of  atmospheric  sub¬ 
micron  particles,  a  systematic  ENA  study  tinder  controlled 
laboratory  conditions  of  high  humidity  (AM  •  90%  to  99.9%)  could 
reveal  growth  functions  for  particular  aerosol  ensembles  and 
elucidate  the  various  stages  and  time  seales  that  H20  molecules 
undergo  in  moist  air  until  they  become  barel/  visible  sa  clusters 
of  >  10*  molecules. 
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01-D4  02-MS  doublets  (Table  1) 

Ewa  excess  water  vapor  absorption 

El 8  computer  ehort-form  for  10*8 

EMF  extreme  high  frequency  range  (30  to  300  GHz) 

HIT RAN  high  resolution  tran- mitilon  model  developed  by  AFGl 

H-bond  molecular  hydrogen  bond 

ID  quanttss  number  identification 

IK  imaginary  part  of 

02'MS  oxygen  microwave  speetrix* 

A£  real  part  of 

WL-*7  atmospheric  transmission  windows 
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ai-a5 

A 

b 

bi-bj 

B(t) 

c 

C 

d 

D(V) 

a 

E 

f 

F* 

F* 

9(f) 

h 

i 

I 

3 

k(t) 

t 

L 


"o 

■Mv) 

M*(V) 

» 

6p 

Q 

r 

RK 


Si,! 

ds 

s 

sa 


? 

e 

n 

x 


coefficient  for  fitting  •  o.  Tabel  5 

Oj-MS  line  coefficients,  Eqs.  (21)  to  (23)  and  Table  L 

total  attenuation,  Eq.  (6) ,  da 

coefficient  for  fitting  data  «  0  »  Table  5 

H2O  line  coefficients,  Eqa.  (24)  and  (25),  Table  2 

second  virlal  coefficient,  Eq.  (37)  *3 

speed  of  light,  Eq.  (3),  km/s 

continue  magnitude,  Eq,  (S3) 

molecular  spacing.  Section  III,  m 

dispersion  spectrum,  Eq.  (13),  ppm 

water  vapor  pressure,  Eq.  (11),  krk 

field  strength,  Eq.  (2) ,  V/m 

relative  hwidity  RH,  Eq.  (45) 

line  ahape  for  D(v) ,  Eq.  (IS),  GHz*1 

line  shape  for  N"(v),  Eq.  (16),  Qtt“^ 

aerosol  mass  growth  factor,  Eq.  (47) 

altitude.  Eg.  (32),  km 

running  spectral  line  index.  Tables  1  and  2 
interference  coefficient*  Eq.  (23) 

-  Eqs.  (1)  and  (3) 

dimer  equilibrium  factor,  Eq.  (39),  s"J 
mean  free  path  length.  Section  IXX,  um 
radio  path  length,  Eq.  (2),  km 
molecular  mass ,  Eq.  (20) 

molecular  number  density,  Eq.  (36),  molec/m3 

refractive  index.  Fig.  13 

real  part  of  n,  Bq.  (SB) 

imaginary  part  of  n.  Fig.  13 

complex  refractivity,  Eq.  (1),  Pf* 

frequency- independent  part  of  N,  Eq.  (12),  ppm 

refraction  spectnmu  Eq.  (27) ,  pfm 

absorption  spectri*.  Eq.  (141 ,  pfxi 

dry  air  pressure,  Eq.  (12),  kFa 

fluctuations  of  p.  Eq.  (34),  kFa 

resonator  quality  factor.  Table  3 

aerosol  particle  radius.  Eq.  (35).  um 

relative  humidity,  Eq.  (45),  % 

ray  path  length,  Eq.  (6) ,  km 

ray  path  coordinate,  Eqs.  (6)  and  (B) 

path  length  increment,  Eq.  (6) ,  km 

line  strength,  Eqs.  (21)  and  (24),  kHz 

aerosol  surface  area,  Eq.  (46),  cm2 

300/T,  normalized  imrera*  T,  Eq.  (11) 

temperature,  Eq.  (11) ,  A 

brightness  temperature,  Eq.  (10) ,  A 

visibility,  Eq.  (52),  km 

air  volume  containing  water  vapor.  Section  IV,  as3 

total  precipi table  water  vapor,  Eq.  (33) .  cm 

total  preeipltable  liquid  water,  Eq.  (49) ,  am 

concentration  exponent  (0),  Eq.  (S3) 

temperature  exponent  (t) ,  Eq.  (S3) 

frequency  exponent  (v) ,  Eq.  (S3) 

attenuation  rata,  Eq.  (4) ,  dB/km 

line  width.  Eg.  (22),  GHz 

Doppler  line  width.  Eg.  (20) ,  kHz 

propagation  constant,  Eq.  (3) 

deviation  from  ideal  9as  behavior,  Eq.  (38) 

propagation  delay  time,  Eq.  (9) ,  ns 

time  constant  for  adsorption.  Section  IV. B.,  s 

time  constant  for  desorption.  Section  IV, B. ,  a 

slant  path  angle.  Section  II,  degree 

optical  wavelength.  Section  III,  n 

frequency,  Eq.  (3) ,  GHz 

line  resonance  frequency,  Eq.  (IS) ,  GHz 

water  vapor  concentration,  Eq.  (36),  9/m3 

path-averaged  value  of  0,  Eq.  (52) ,  9/m3 

fluctuations  of  p,  Bq.  (34),  9/a3 

transmittance ,  Eq.  (7) 

phase  rate,  Eq.  (S),  rad/km 

total  phase  delay,  Eq.  (8),  rad 

sum  of,  Eqs.  (13)  and  (14) 


Subscripts 

A  aerosol 

D  dimer 

d  dry 

(  absorption  line 

L  Lo'  ♦"tzian 


% 


I* 


*  H 


Acronyms 


Ara 

at 

cat 


Air  Force  Geophysical  Laboratory 
condensation  nuclei 
clood  condensation  nuclei 


maximum 

integer,  cluster  sin 
Initial  valua 
saturation 
water  vapor 
liquid  watar 
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S  EWA 

U  Optical  wavelength  ( i  -  0.5$  ue) 
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ABSORPTION  OF  MILLIMETER  TO  SUBMILUMETER 
WAVES  BY  ATMOSPHERIC  WATER  MOLECL'LES 


Dtprnrtmtm  of  Pkpaa 
Umvomiy  of  Coferado 
BtmUtr.  Co*m*H3» 


Using  the  values  of  the  rotational  Molecular  parameters 
(including  centrifugal  distortion  term)  of  the  Bj,40 
oolecula.  which  c an  explain  12  observed  transitions  below 
800  GHs,  all  rotational  energy  levels  with  normalised 
Boltzmann  factors  larger  than  5x10"*  at  300*8  are  calcu¬ 
lated.  Probabilities  of  all  possible  elsctrlc  dipole 
transitions  among  these  states.  2277  lines ,  era  calculated 
using  the  a  1  gen  functions  thus  obtained,  and  the  permanent 
electric  dipole  moment  of  1.8546  Debye.  AssiaUng  the 
single  and  full  Lorentsian  line  forma,  we  calculated  the 
absorption  coefficient  for  millimeter  to  sub mil lime*  - 
region.  Our  result,  using  the  single  term  Lor eotr Lao  line 
form,  agrees  quite  well  with  experiment  for  1  Torr  of 
water  vapor  in  760  Torr  air  at  300*8. 

Key  words:  water  vapor,  allllmeter  waves,  submill lmeter 
waves,  excess  absorption,  Lorentsian  line  form,  wlodow 

regions. 


Introduce  loti 

Burch  (1)  observed  absorption  of  millimeter  to  sub- 
millimeter  waves  oy  Mist  air.  particularly  l«  the  “window 
regions",  ind  found  that  the  experimental  values  were 
larger  than  the  values  calculated  lor  those  resonance  lines 
nelov  1  THx  bv  a  factor  of  about  two.  Such  eacess  attenua¬ 
tion  was  confirmed  by  both  laboratory  end  field  experiments 
(2)  and  has  been  recognised  in  mmerlcal  calculations  as  an 
empirical  correction  term  (3) 

«u  -  *.Z3x10"V<300/T)J-1<v/30)J  «/■» 

where  p  la  the.alr  pressure  in  kPa,  o  la  tha  water  vapor 
danalcy  la  g/a  ,  T  la  tha  temperatura  in  *8,  and  v  la  tha 
frequency  in  GHx. 

McClatchy  and  others  (4)  gave  a  list  of  1838  B?0  lines 
balow  31  THx,  baaed  on  a  rigid  body  aodel  of  tha  molacula, 
and  an  seal- empirical  formula  for  the  line  width: 

4  -  bj(p  +  *.a0pw)(300/T)0'6 

where  p  is  the  partial  preaaora  of  tho  water  vapor.  Llaba 
(5)  svsmserlsed  tha  ealculatod  excess  absorption  duo  to 
higher  frequency  lines  as 

.  l.S*10'5pe<»0/T)(v/»)2 

Crane  (6)  reviewed  the  situation  presenting  more 
r scent  experimental  dote. 


The  theory  of  rotational  states  of  aeywetrlc  top 
molecules  is  well  established  (7,6).  For  the  water 
molecule  Hj^o,  and  Its  isotope  aolecules,  DeLucla  and 
others  (9.10)  found  seta  of  waives  for  11  molecular 
potMeters,  including  three  rotational  cons  tenet  and  alghe 
centrifugal  distortion  constants,  by  which  all  transition 
frequencies  below  800  GHs  they  observed,  could  be 
explained*  We  singly  took  their  values  of  the  nolecular 
parameter*,  calculated  energy  levels  up  to  J  •  16,  and 
found  chat  all  states  with  normalised  Boltxnarm  factors 
larger  than  5x10*®  at  300*8  (8x10*®  at  313*8)  wera  counted 
in.  Ihe  hfe  splitting*  wera  not  calculated,  but  the 
statistical  weight*  due  to  the  proton  spins  wore  taken 
into  account.  Hit  frequencies  and  intensities  of  all  2277 
electric  dipole  transitions  amont  these  rotational  states 
were  calculated  using  cha  eigenfunctions  obtained  above 
and  tha  accepted  vslue  of  the  dipole  moment,  1.8546  Dehve. 

Ua  then  assumed  that  each  line  had  the  Lorentsian 
form,  and  added  them  up  to  obtain  the  absorption 
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where  is  cha  nusber  of  aolecules  in  the  s-state  per  unit 
volume,  u  la  a  apace  component  of  the  electric  dipole 
moment,  and  v  la  tha  transition  frequency  between  tha  e- 

aad  b-etatee. 

The  Lorentsian  form  foccor  given  in  the  above  formula 
la  called  tha  full  Lorentsian  (abbreviated  as  f-L).  Since 
the  second  term,  due  to  the  induced  emission  of  photons, 
oay  be  suppressed  by  scatterings.  It  is  worthwhile  to  try 
tha  single  Lorentsian  (abbreviated  as  S-L)  form  factor  in 
which  tha  ascend  term  is  neglect  ad.  We  triad  both  F-L  and 
S-L  form  factors. 

We  asaiaad  that  tha  transit  ion  frequencies  are  all 
independent  of  praacure  and  the  width  parameter  la 


»  *  jig*  <*  ♦  *.*0|>w)(»om0-s  au 


as  observed  experimentally  for  soma  lines.  No  quanta 
number  dependence  was  assumed. 


Result 

Our  results  for  300*8,  p  •  760  Torr,  and  Pw  •  1  Torr 
are  shown  in  Table  1  for  some  frequencies  up  to  1  THx ,  and 
compared  with  the  aodel  attenuation  formula  (5)  which  fits 
available  experimental  data.  We  sea  chat  F-L  is  b*L<nr 
experimental  values ,  but  S-L  result  agrees  quite  well  with 
experiment . 

As  la  known  alraady,  the  26  lines  below  1  Tils 
contribute  only  falf  of  the  absorption  coefficient  in  that 
region,  out  the  remalninx  nail,  tne  so-called  exceen 
absorption,  is  simply  tne  accumulated  absorptions  ot  the 
higher  frequence  transitions  in  the  nrrsent  theory  with  the 
s-L  form  lector.  Figure  1  shown  our  lomnuter  plotting  ot 
rhe  water  spectrum  below  2  THx  using  tne  S-L  fora  factor, 
ind  it  fa  almost  identical  to  tha  seal -empirical  ancctrua, 
except  for  the  oxvgen  aosorpcions  which  are  not  considered 
here . 

We  also  calculated  the  absorption  coefficient  at 
300*8,  p«  760  Torr,  and  py  -  10  Torr  with  the  S-L  form 
factor  and  saw  that 

log  *  HO  Torr)  -  luge  (l  Torr)  ♦  1.0227 

hoided  for  V>  (.Hi  through  1  THx,  except  near  resonant 
regions. 

The  tecneratura  dependence  appears  through  both 
Boltzmann  factors  and  tne  width  parameter.  Calculation, 
with  the  S-L  form  factor,  was  sade  for  273*8,  p  -  760  Torr, 
and  -  l  Torr,  and  we  touna  chat 

log  a  (273*8)  -  logo  (300*8)  *  0.103 

hoided  for  the  low  frequency  region  up  to  150  GHs,  and 
that  the  additional  term  went  up  to  0.108  around  200  GHs. 

Including  the  absorptions  by  tha  atmospheric  oxygen 
molecules  calculations  were  made  for  temperatures  up  to 
1000*8,  and  wa  found  that  tha  absorption  in  tha  200  GHs 
region  would  be  quite  low. 
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This  Interaction  produces  the  magnetic  dipole  transition* 
between  rotational  state*  given  by  (2).  The  matrix  elements 
of  (3),  which  give  intensities  of  the  transitions,  are 
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In  addition  to  the  water  molecule  which  we  considered  in  a 
previous  paper,  the  oxygen  molecule  la  taken  Into  account 
to  express  theoretically  the  absorption  coefficients  of 
micro-,  millimeter- ,  and  submlllimecer-vaves  by  air  at 
given  temperature  and  preaaurc.  The  theory  la  applied  to 
each  layer  of  the  U  S  Standard  Atmosphere,  and  the  result- 
ing  absorption  coefficients  are  Integrated  for  altitude 
from  0,  2,  and  4  km  to  80  km  to  obtain  the  transparency  of 
earth's  atmosphere. 
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Key  words:  transparency,  absorption  coefficient,  oxygen 
molecule,  water  molecule,  U  S  Standard  Atmosphere. 
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In  a  previous  paper  (l)  we  calculated  all  rotational 
states  up  to  J«  lb  of  the  water  molecule,  using  known  rota¬ 
tional  and  centrifugal  distortion  constants,  and  all  possi¬ 
ble  electric  dioolc  transit  Ions  among  these  states,  2277 
line:..  Tt  was  shown  that  the  ah»orntion  of  micro-,  milli¬ 
meter,  and  suDmillimeter-w.-ivrs  by  1  Torr  of  water  vapor  in 
air  of  760  Torr,  300  K  could  be  explained  if  we  assumed 
that  each  of  theaa  lines  had  Loren trian  line  shape  with  a 
common  width  parametar  determined  by  experiment  on  some  of 
them.  It  was  saen  that  th#  single-Loren tizlan  line  shape 
gave  almost  complete  agreement  with  empirical  values  of  the 
absorption  coefficients,  while  Che  full-Lorentxiam  line 
shape  gave  calculated  values  about  1/3  less  than  the  empir¬ 
ical  values,  which  was  actually  not  so  bad  conaldarlng 
large  uncertainties  in  the  empirical  values. 

In  the  present  calculation  we  take  the  full-Lorenttlen 
line  shape  since  It  has  better  theoretical  foundation  chan 
the  alngle-Lorentzlaa  line  shape. 

Rotational  Transitions  of  Che  Oxygen  Molecule  (X 

The  microwave  spectrum  of  the  oxygen  molecule,  , 
in  its  electronic  ground  state,  X  ^IJ,  has  been  measured 
carefully  (2-4),  and  th#  most  recent ^measurement  by  Ends 
and  Hiroshima  (5)  is  shown  to  be  reproduced  by  the  follow¬ 
ing  Hamiltonian 

H  .  BR2  -  DN4  +  <2/3)<»<3S,2-S2)  ♦  (1/2)»0(ICJ(3SC2- 
S2)+(3S|.2-S2)N2)  ♦  (1/2)»od((I4(3S?2-S2)+(3Sc2-S2) 
N4))  +  yR-S  ♦  YdK2R-S  (1) 

with  propter  values  for  cha  parameters,  8,  D,  1,  X  ,  1  ,  y, 

end  Tjj. 

The  resulting  eigenstates  can  be  expressed  in  terms 
of  the  eigenstates  of  *2,  s2,  J2,  and  Jg  operators,  {NSJM>, 
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-here  each  eigneatate  is  expressed  as  tnJHx.  The  value  of 
the  S  quantum  number  is  Uxcd  J*  1. 

The  wldthe  of  the  microwave  tines  of  the  oxvgen  mole¬ 
cule  are  measured  (7)  to  be  2. j  MHs/Torr  and  nearly  indepen¬ 
dent  of  the  quantum  number*.  The  widths  ot  the  430  and  760 
^Hx  lines  are  not  measured  vet.  hut  we  simplv  assume  that 
ait  oxvgen  lines  have  the  tull-Lorentxian  line  shape  vitn 
comon  width  2.3»'300/T  MH*/Torr.  The  absorption  coefficient 
for  a  typical  case  of  760  Torr,  300  l  alt  with  l  Torr  of 
water  vapor  is  calculated  considering  the  oxygen  and  water 
molecules  in  this  way.  The  result  is  shown  ir  Fig.  1. 


Figure  1.  Computer  plotting  of  the  absorption 
coefficiant  by  air  of  760  Torr,  300  K  with 
1  Torr  water  vapor* 


Transparency  of  the  U  S  Standard  Atmosphere 


where  S  is  not  shown  explicitly  since  It  is  fixed  as  l,  and 
a,  and  are  roof f lclrncs  determined  for  each  J. 

rite  interaction  ot  tnc  magnetic  moawwt  of  this  mole¬ 
cule  with  external  field  B  is  given  (6)  by 

"  *  -  .(it. *.».*  t:  is  J  «  «  )  -  V»-»)  (3) 

'.»•  n  *■  l»  -V  -  -  V  J* 

-here  u ^  is  the  bohr  magneton,  <.  is  the  molecular  axis, 

Sj  -  Sf  *  iSa.  and 


Typical  values  of  temperature  and  pressure  at  each 
altitude  of  earth's  atmosphere  are  given  in  the  U  S  Stand¬ 
ard  Atmosphere.  We  apply  the  above  explained  method  of 
computing  tnc  absorption  cowtficicnt  to  each  I  aver  ot  the 
If.  S.  Standard  Atmosphere  and  then  integrated  them  over  the 
altitude  up  to  80  km  (top  of  the  mesosphere)  to  obtain  the 
transparency, 

°0  km 

transparency  •  expt-  i  Jh)  (6) 

x 

The  lower  limit  of  the  integral,  x,  la  chosen  to  be  0,  2, 
and  4  km  in  the  computations,  results  of  which  are  shown  in 
Fig.  2a,  2b,  and  2c,  respectively. 
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Figure  2a.  Transparency  of  Che  D  S  Standard 
A taw* ph era  it  0  ka  altitude  la  0  to  1  THr . 


Figure  2b.  Transparency  of  the  U  S  Standard 
Atmosphere  at  2  km  altitude  in  0  to  1  THz. 


. . . 


Figure  2c.  Transparency  of  the  U  S  Standard 
Atmosphere  at  4  km  altitude  in  0  to  1  THz. 
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It  ii  luuiid,  In  this  calculation*  that  earth's 
atmosphere  contains  20 .951  of  1402  sC  »11  altitudes.  The 
saturation  pressure  of  water  vapor  dacraases  rapidly  as  the 
altitude  increases  as  shown  In  Tabls  1.  Since  it  la  so  snail 


Table  I.  Saturation  Vapor  Pressure  (Tort)  as 
a  Function  of  Altitude  in  U  S  Standard  Atmosphere. 


ka  012343*7 

Torr  12.79  8.33  3.29  3.14  1.91  1.09  0.58  0.29 


K  ¥  10  11-24  23-27  2A-2H  ,'i 

0.13  0 .08  0.03  0.02  0.03  0.04  u.OS 


at  thr  '!  r«»popaujic(20  im)  w.*  assume  that  wat*-r  v.npor  exist 
only  in  altitudes  below  20  km.  and  chat  the  hu-iditv  *tavs 
constant  iron  0  to  20  km,  to  limplify  conputatlon.  The 
huniditv  In  this  sense  is  taken  as  02,  10" .  302,  and  1007 
in  Figs.  2.i  through  2c. 
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Absolute  attenuation  rates  due  to  water  vapor  and  moist 
nitrogen  have  been  Measured  In  tbe  laboratory  at  138  CHs, 
282  and  300  K,  pressures  up  to  1.3  atm,  and  relative  humid¬ 
ities  fro*  60  to  100  percent.  The  coaepucer-cont rolled 
eessuring  system  is  comprised  of  a  millimeter  wave  reson¬ 
ance  spectrometer  (0.15  km  effective  path  length)  and  a 
humidity  simulator.  Several  shortcomings  of  earlier  meas¬ 
urement  attempts  have  been  rectified.  The  data  are  inter¬ 
preted  as  a  water  vapor  continuum  spectrum  consisting  of 
two  terms,  namely  strong  self-broadening  (H3O»H20)  plus 
fore lgn-g as-broadening  (H^O-Nj)  contributions.  Implica¬ 
tions  of  the  new  results  for  modeling  atmospheric  DIF 
window  transparencies  and  for  revising  established  H20  line 
broadening  theory  are  discussed. 

Key  Words:  Atmospheric  radio  wave  propagation;  high- 

humidity  spectrometer;  millimeter  wavelength 
region;  water  vapor  continuum  absorption 


•The  work  was  supported  In  part  by  tbe  U.S.  Army  Research 
Office  under  Contract  AR0  101-83. 


(dimers,  clusters,  etc.)  might  be  involved.  Results  fora, 
froe  controlled  laboratory  experiments  in  chs  140  GHz 
window  are  reported  which  help  to  clarify  the  watcT  vapor 
excess  problem. 

Water  Vapor  Millimeter  Wave  Spectroscopy 

Water  in  both  vepor  and  liquid  states  obscures  the 
millimeter  wave  transparency  of  the  troposphere.  Frequency 
dependence  and  magnitude  of  H20  attenuation  are  both 
distinctly  different  for  the  same  (1  am)  absorber  thickness 
in  wstsr,  rain,  suspended  droplet,  moist  sir,  and  pure 
water  vapor  form,  as  Illustrated  In  Flgyce  l  (2). 


FREQUENCY  f,  GHz 

Figure  l.  Normalized  specific  H2O  attenuation  over 
the  frequency  range  5  to  1000  GHz  for  five  phase  states: 
water  Q>  rain  QX  suspended  hydrometeors  QX  moist 
sir  0^  (numbers  shown  indicate  selected  H20  attenua¬ 
tion  peaks),  and  pure  water  vepor  QX 
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Introduction 

Millimeter  weve  propagation  through  the  atmosphere  is 
affected  by  meteorological  phenomena  such  as  humidity,  fog, 
clouds,  sad  precipitation.  System  applications  require  a 
detailed  knowledge  of  atmospheric  propagation  limit at Iona 
and  their  dependency  upon  relevant  meteorological  parame¬ 
ters  .  The  Influence  of  tbe  Intervening  atmosphere  upon 
radio  waves  can  be  quantified  by  s  propagation  modal.  In 
daar  sir,  water  vapor  and  oxygen  absorption  causa  fre¬ 
quency-dependant  signal  attenuation,  propagation  delay,  ray 
bending,  and  mad  lisa  noise.  Prediction  of  these  propagation 
affects  la  accomplished  by  s  broadband  modal  of  the  complex 
refract ivlty  N  (1).  To  cover  frequencies  up  to  300  CHx,  a 
reliable  expression  for  N  requires  more  than  200  coeffi¬ 
cients.  Tbe  refractlvity  N  is  formulated  la  units  of  parts 
per  million  (ppm)  sad  consists  of  throo  components 

»  -  I!©  ♦  *><!)  ♦  JH"(D  PPm 

where  W  is  the  frequency  independent  refractlvity,  and 
0(f)  and  !i**(f )  are  frequency-dependent  dispersion  and 
absorption  terms  arising  from  the  molecular  spectra  of  H20 
end  0j. 

The  specific  rates  of  tbv  propagation  eedlia  are  power 
attenuation 

o  •  0.1820flT<f)  dB/km  (1) 

end  daisy 

6  -  3.336(11  ♦  0(f)]  pa/km  ,  (2) 

o 

where  the  frequency  f  la  la  units  of  GHs. 

This  paper  addresses  ailllmeter-vave  water  vapor  at¬ 
tenuation  rates.  There  has  been  a  long-standing  discrep¬ 
ancy  between  theoretical  spectroscopic  predictions  and 
experimental  data  obtained  from  both  laboratory  and  field 
measurements  (3,4).  Marked  differences  occur  In  the  window 
ranges  located  around  35,  90,  140,  and  220  CDs  (2),  where 
experimental  attenuation  exceeds  the  theoretical  valuaa. 
Since  excess  factors  between  2  end  3  ere  involved,  it  has 
been  speculated  Chet  unidentified  absorbers  related  to  H20 


Water  vapor  attenuation  depends  on  absolute  tumidity 

»  -  7.219*9  ,/»5  (3) 

where  e  Is  the  partial  vapor  pressure  in  Ufa  -  10  mb  and 
e  -  300/T  is  a  relative  Inverse  tasperature  parameter. 
Water  vapor  pressure  e  Is  limited  by  relative  tumidity 

U  -  *l.S1.9_5lOt,-M*e-10>  <  1001  .  (*) 

The  attenuation  has  Its  origin  in  absorption  Unas  of 
tbe  Bolecule  H20.  la  principle,  a*  can  be  evaluated  from 
a  llae-by-llne  ■unmet ion  of  all  line  features,  each  con¬ 
tribution  at  a  given  frequency  as  prescribed  by  a  line 
shape  function  F(f).  Seven  reported  H20  line  data  bases 
(LB)  are  indexed  in  Table  I. 


Table  I.  Available  date  baaee  for  1)0  absorption  limns 


10 

Frequency 

cutoff 

Number 
of  llmeo  (n) 

Remark 

Ref. 

LSI 

ns 

337 

38350 

A7CL  Tape 

3 

U 2 

126 

17201 

French  version 

6 

LB3 

30 

2277 

locational  Spacer. 

X  7 

L84 

3 

LSI 

Radio  Astronomy 

8 

105 

1 

30(*O 

RHF  approximation 

1,23 

104 

1.1 

17  c 

Russian  version 

18 

107 

23 

535 

Rotational  spacer. 

XX  27 

The  extreme  low  frequency  vtnge  of  about  180  stronger  (l.e., 
at  vQ  more  than  10s  times  the  value  of  a*  at  100  Ola)  lines 
below  12  Til  make  contributions  In  the  CRT  range.  Six  dif¬ 
ferent  standard  line  shape  functions  F(f)  are  available  to 
assess  these  contribution!.  Typically,  the  shape  function! 
are  governed  by  two  spectroscopic  parameters:  line  center 
freeueney  vQ  and  width  y.  An  additional  parameter  «.  was 
Introduced  for  covering  a  wide  spectral  range,  which  u 
defined  by 

vB  •  2kT/he  -  12310/6  OU  (3) 

(4b,  4h).  When  these  shape  functions  are  normalised  ta 
unity  at  the  center  vQV  the  functional  dependence*  given 
in  Table  II  are  obtained. 
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Tabic  II,  Normalized  pressure-broadened  line  shape 
functions 


Table  111.  Theoretical  EHF  water  vapor  continuum 
coefficients  (Eq.6) 


Normalization:  for  a  single  line. 


ID* 

Fix) 

«»io3 

(Ex.  7) 

SL 

•2«  (»2+(l-l)2]'1 

1.9 

Ft 

»2x  (I.2.!!-.)2)-1  -t»2+u-.)2l_1) 

-0.001 

CR 

»2*2  uV-Hl-.2)2/*!*1  »  (Ft) 

-0.001 

W 

•2x2tfx2+a-x)2r1  ♦r«2+a>*)2rl) 

637 

MCR 

(C*)z 

8.2 

MTU 

(VV)  1 1 xnh(bx) /tab  b)/z 

X6.6 

SL  •  shore  Lorentzian,  Ft  -  full  Lo  rents  las , 

GR  -  Cross,  VW  -  Van  Vleck-Weisskopf , 

MGR  -  modified  CR,  MW  -  modified  W  (10). 

The  intensity  distribution  of  a  strong  HiO  line  vas 
traced  with  the  shape  functions  listed  in  Table  11  and 
the  example  la  exhibited  in  Figure  2.  Within  the  line  core, 
that  is  F(z)  >_  10'“  or  t  2  0.8  to  1.2,  ell  shapes  ere  alike 
and  it  Is  Inconsequential  which  one  Is  used.  On  the  other 
hand,  far- wing  intensities  depend  very  much  on  the  chosen 
frequency  distribution  function.  For  exompls,  the  relativs 
contribution  of  the  line  shown  in  Figure  2  can  vary  at 
100  CHs  over  a  range  l(HCR) :  75(51.) . 


The  functions  were  tested  for  finite  integrated 
(f  •  0  to  •)  line  absorption  to  decide  which  shape  is  most 
plauslbls.  A  limited-range  maerlcal  Integration  resulted 
in  residuals, 

1000 

6  •  (l/wa)/r(t)d*  -  1,  (7) 

0.001 

and  It  follows  (son  Tab  la  IX)  that  some  deviate  consider¬ 
ably  from  the  ideal  6-0.  Based  on  Eq.(7),  the  shapes  W, 
MW  and  MCR  can  not  be  considered  valid  intensity  descrip¬ 
tor  of  e  preeeure-broedened  line  over  e  wide  frequency 
range  (22) — and  the  othera  are  probably  overt axad  by 
•^Rtctlng  thorn  to  predict  eight  significant  orders  of 
sttsnuatlon  rate. 

Water  Vapor  Continuum  Absorption 

The  attenuation  rate  a  due  to  a  single  line  is  calcu¬ 
lated  with  Eq.(l),  where  the  product  of  strength  5  sod 
shape  F  defines 

»T(f)  -  SF(f )  ppm  .  (8) 

The  line  properties  ere  formulated  In  Ref.(l)  at  0  -  1 
(300  K)  by  a  set  of  spectroscopic  coefficients  bi  to  b-. 

In  the  low  frequency  far  wing  (y«f«v  ),  any  of  the  theo¬ 
retical  line  shapes  In  Table  II  is  proportional  to  the 
width  y;  however.  Individual  continuum  shape  approximations 
F  (f  *  0)  - 
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Figure  2.  Nornellz.d  preeeure-broedened  line  xhape* 
F(z)  (eee  Table  II)  applied  to  a  single  H20  line. 

The  line  is  actually  centered  at  v0  •  1669.907  GHz 
with  a  peak  attenuation  rate  a0  -  2 . 3  «  10s  dB/taa 
In  sea  level  air  <9  «  1,  v  •  10  g/a3). 


*fy/v’  (Gl.FL).  T/v‘  [  SL  1. 

predict  succu.lv.ly  hlghm  contrlbucloa.  («M  Fl«.  2, 
Tabl.  III). 

la  octe  to  fouulat.  ch.  conclu.  .  Of.  6.  Tabl.  ill) 
u  a  m  e{  u  li»  mu  ducxlb.A  by  altbu  Tina 
.Upu,  ch.  following  .uuptlmi  hm  to  b.  ud.  .but 
.crutch  u4  wldcb  of  u  "iqulnlat'  Uu  mbamm  1  TU: 

s  -  b,u3-s  up  tbjtl-e))  •  bju1'5  uu  (io) 

where  bj  *  2.0,  and 

y  ■  bg(m*e  ♦  p)0  CBs  (U) 

Egs.(l),  (I)  to  (IX)  Xeed  to  the  express** 

•  •  0.7312 •(...  ♦  pie1-5  .  I  OibVv5).  .  02) 

*  t  u  1  "  4 


The  continuum  spectrum  consists  of  i  self  kmdta 
(■  a2)  and  a  foreign-gas '  broadening  term  (-  ep) , 


•e  ■  (hga*  ♦  kfep)e^  dB/taa 

GoeffielaBts  and  temperature  aiponent  ere  i 
shape  theory  os  follows: 

k.  -  O.naf2  l  »lbl/vj)t  -  »kf  . 


The  EHF  (30  to  30 0  GHz)  propagation  model  developed 
by  us  (1)  uses  for  reasons  of  computational  economy  30 
local  H?0  lines  below  1  TRx  (LBS,  Table  I).  A  theoretical 
assessment  of  EHF  attenuation  +  a  wee  accomplished 

on  the  bos la  of  LB3  by  subtracting  the  line  contributions 
below  l  THx  from  all  2277.  The  far-ving  contributions 
frM  pres sure- broadened  lines  above  1  Tllz  are  then  fitted 
below  300  CHs  to  a  water  vapor  continuum  absorption 

cc  -  k(p/101)vV(f/100)U  d«/U.  (6) 

where  p  is  dry  air  pressure  in  kPa.  One  recognizes  ex¬ 
ponents  for  vapor  concentration  (x),  temperature  (y),  and 
frequency  (u)  dependences.  Applying  various  shape  func¬ 
tions  (Table  XX)  to  LB3  led  to  the  results  summer iced  in 
Table  XIX.  As  expected,  the  scale  factor  k  varies  over  a 
wide  range.  Also,  one  notices  that  all  three  exponents 
ere  influenced  by  the  selected  shape  function. 


y<s)  -  2.3 
The  reels 

■■‘A 


T'(»>  -  1.5 


C.SO(Air) , 


is  the  pressure-broadening  efficiency  of  pure  water  vapor 
with  referenca  to  dry  air  or  nitrogen*  as  determined  from 
measurements  in  line  cores  (IS). 

There  Is  ample  evidence  (6a,  10,  12)  that  in  the  EHF 
range,  besides  local  lines,  the  low  frequency  wing  of  the 
complete  self -broadened  rotational  B|0  spectrum  Is  much 
stronger  chan  predicted  by  Eg.  (12).  1  careful  review  of 

HjO  continuum  absorption  by  Burch  (10) •  Imeludlng  all  of 
his  reported  date,  uncovered  increases  of  tha  broadening 
efficiency  In  the  far  wings  of  lines  ranging  between  a  • 
20  at  J00  CH*  end  ■  ■  SO  at  30  CN>  (see  Later,  Table  IV). 


-  .  *  ,**  -*» 


Tht  laboratory  experiment  described  In  tha  next  sac  t ion 
was  prepared  to  measure  absolute  atcanuaeion  ratas  o  at 
thraa  EHF  window  frequencies  (1.*.,  92,  138,  and  220  GHz ) 
under  conditions  ranching  watar  vapor  saturation.  Ab¬ 
sorption  by  nalghborlng  liras  is  snail,  so  that  to  first 
ordar  a  «  a  .  Attanuation  data  ara  expected  to  behave  as 
pr  ad  let  ad  by  Eq.(13).  Tha  ongoing  experimental  studlas  ara 
alaad  at  sacurlng  accurata  coafflciaats  k  ,kf,  and  y ,  as 
wall  as  looking  for  possible  anomalies  not  supported  by  tha 
pressure- broadened  rotational  spactnas  of  tha  HjO  aolecule. 
Tha  last  s act ion  of  this  papar  praaanca  first  attanuation 
raaulta  obtained  at  138  GHz  and  discusaaa  thair  implica- 
tloos  to  EHF  propagation  nodaling. 

Millimeter  Wave  Spectrometer 


Controlled  experiments  that  simulate  atmospheric 
conditions  provide  cast  caaaa  for  studying  specific  contri¬ 
butions  to  If  in  isolation.  Aaanaaesr.ca  of  basic  physlenl 
principles  uadsrlylng  the  attenuation  rats  a  ara  difficult 
to  make  free  eaaauraeants  la  tha  actual  ateoaphara.  la* 
liability,  praclsion,  and  seals  of  supporting  metaorologl- 
cal  dots  cn^-nel  ss  tha  quality  of  eoat  f laid  observations. 
Tha  expected  EHF  window  attanuation  la  low  la  a  spectro¬ 
scopic  sanaa  (<10  dB/ke)  and  raquiraa  long  path  lengths 
(>100  a)  for  s  reliable  detect  Ion.  A  saal-confocal  Fabry- 
Pdrot  resonator  with  a  high  Q  value  (137,000  at  138  GHz) 
offered  a  good  coeproaiso  between  compactness  and  dataction 
sensitivity.  Tha  cneputar-controllad  aaaauring  system 
consists  of  the  ail  lime  tar  wave  resonance  apace  remoter  and 
a  humidity  simulator.  An  Insulated  box  contains  s  high- 
vacuum  stainless  steal  vassal  which  houses  s  temperature- 
controlled  mini- lake  (10  cm  across)  end  tbs  radio  test  path 
(see  Figure  3).  Signals  sc  frequencies  betvtai  30  sad  300 
GHz  can  be  employed  over  equivalent  f relapses  path  lengths 
between  500  sad  50  a;  temperatures  era  controlled  to  batter 
than  1M0O  of  a  degree  Celsius;  pressure  ranges  over  seven 
orders  of  magnitude  (102  to  10~s  kPa);  and  relative  humid¬ 
ity  is  varied  between  0  sod  100  percent. 

Schcmatical  diagrams  of  physical  (Figure  3)  and  elec¬ 
tronic  (Figure  4)  arrangements  comvsy  an  overview  of  Che 
experiment.  Quasi- static  diffusion  nixing  was  chosen  to 
generate  well  defined  tumidities  inside  the  spectrometer 
cell.  A  temperature  controlled  water  reservoir  servos  ae 
the  vapor  source.  Elect ropolished  stainless  steel  was  used 
exclusively  as  construction  material.  Various  hydrophobic 
coatings  were  studied  as  possible  means  for  neutralizing 
the  absorption/ desorption  cycle  on  surfaces  exposed  to  water 
vapor  <4d),  but  were  abandoned  in  favor  of  slightly  heating 
the  mirrors.  Four  fast-responding  (r  <  Is)  temperature 
sensors  Inside  the  ceil  diagnose  any  disturbance  of  the  gas 
equilibrium.  Data  acquisition  is  computer-control  led. 
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Figure  3.  Cross-sectlcm  of  millimeter 
vovo  resonance  spectrometer  and  built-in 
humidity  generator.  Note  three  different, 
extremely  stable  temperature  sones,  Tj (0j ) 
to  TjCOy).  Relative  humidity  inside  the 
cell  is  (Eq.  4):  log  (RH)  - 

2-9.834  (0,  -  02)  «•*  5  (log©!  -  log©  )  X. 
Absolute  faaldity  v  is  shown  in  the  diagram. 


Figure  4.  Schematic  of  high-humidity  millimeter  wave  spec¬ 
trometer  end  computer-controlled  date  acquisition  system. 


User- friendly  software  allows  the  choice  of  several  modes 
of  continuous  operation  over  time  scales  from  one  minute  to 
many  days.  The  resonance  signal  that  carries  the  absorp¬ 
tion  information  is  detected  with  s  digital  waveform  ana¬ 
lyzer.  The  analyzer  performs  the  following  tasks  under 
computer  control:  averaging  (1  to  lOOx),  baseline  normal¬ 
ization,  peak  value  (P-P)  and  bandwidth  (BU)  meaeurwents. 


Dataction  Principles 

The  dimensionless  complex  refractivlty  K  Is  a  macro¬ 
scopic  measure  of  tha  Interaction  between  gee  molecules  and 
waves.  Real  and  imaginary  parte  are  related  by  the 
Kramer s-Kr on ig  formula.  The  response  curve  of  a  gas-filled 
resonator  is 

•  (f)  -  {<fl  +  ♦  H’.LO'Vll}'1.  (12) 

whers  Q  is  the  loaded  Q-value,  f  is  Che  resonance  frequency 
of  the  evacuated  resonator,  and  a  Q-value  of  the  ges  is  de¬ 
fined  by 


Two  modes  of  spectrometer  operation  follow  from  Eq.(l7) 
vhan  the  evacuated  resonator  is  slowly  subjected  to  s  gas 
pressure  P: 

Absorption  mode.  The  resonator  is  swept  periodically 
across  its  resonance  f«  in  order  to  measure  a.  The 
shape  of  the  reeponse  cum,  Eq.(17),  appears  as  a  pulse 
series  a(t)  when  rectified  with  tha  dataction  law  £.  The 
decrease  In  the  maximum  amplitude  a  (F  •  0)  at  to 
«li».  »,(P)  1.  .«tu«d  by  .  path  d«2«ctor.  TS.  itt  (mutton 
rate  is  obtained  from 

.  -  -ll/t,  4*/k»  .  (19) 

HtMT*  ch*  .Uectlr.  poth  length  to  |tna  by 

l,  -  ..771  •  10*5Q/f,  ka  .  (20) 

From  s(t) ,  the  peak  dateccor  records  an  absorption  preeaure 
profile 

•l(P)  -  .^{[.(pn^/t.MJJ  ♦  l}6  ,  (21) 

where  the  total  pressure  P  •  e  ♦  p  la  measured  and  the 
expected  attenuation  rate  a  Is  predicted  by  Eq.(l3).  A 
digital  waveform  processor  (sea  Figure  4)  resolves  amplitude 
changes  to  about  0.1  pcrcant,  allowing  s  detection  sensi¬ 
tivity  of  (16) 

%•«  i 

(..5.,  rR  •  138  Cllz,  0  -  4)7000,  »ata  «  0.05  dB/ka).  (22) 

Refraction  rod..  To  au.ur.  N' ,  th.  thtft  In  rcunnn 
froqu.ncy  by  o  .Mil  .mint  to  f„  u  tracked  and  c.fractton 
Is  given  by 

»'  •  loV,  -  f!J(P)  ]/*K  pp..  (23) 

Frsqusncy  changes  of  f~  can  be  compensated  with  the  tuning 
Micrometer,  allowing  s' resolution  of  N*  0.6  ppm. 
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Specification  for  the  136  GHz  Test  Serlgi 

Resonator.  Fabry-Mrot  rtflmloft-typ*,  ml-conf ocal 
arrangement,  10  cm  mirror  diaaeter,  20  c a  (pacing;  nicro- 
Mtar  timing:  0.3175  am/ cum  with  1.3  «  10  am  resolution; 
loaded  Q  value  at  f  -  138.2  CHx:  q  -  437,000;  effective 
path  length  (Eq.20):  L_  •  0.151  tat;  Freanel  number:  6;  re- 
spooae  tine  constant:  Q/»£_  -  l.Olua;  coupling  coeffi¬ 
cient:  0.053;  coupling  hole:  circular,  0.65  an  disaster, 
0.075  an  double  Mylar  vacuum/ pressure  teal. 

Excitation.  Linear  11  y  frequency  modulated  klyatroo 
with  *  tab  Hired  center  frequency;  (weep  rate:  50  Hr,  sweep 
width  7  MBs  or  19600  us;  reflected  signal  off-resonance: 

150  aV. 

Detection.  Ohbiaeed  Schottky  diode  (8  -  2)  with 
P  recta  ion  attenuator  for  dB-calihratioa;  resonance  signal 
peak  sensitivity  (Eq.22):  t  0.03  dl/te;  accuracy:  better 
then  ±  10  percent  for  »>1  dl/faa. 

Pet  a  tun.  An  example  of  aee  airring  water  vapor  attenua¬ 
tion  By  by  dec  acting  the  resonance  peak  signal  a;  (Pk)  with 
reference  to  vacina  (a  )  la  shown  la  Figure  5.  Relative 
humidity  wee  Increased9  from  RH  -  90  to  100  percent  over  a 
two-hogr  period  by  raising  the  water  toperature  from  6.8 
to  8.4  C.  Condensing  water  vapor  added  several  tenths  of 
a  degree  to  the  cell  temperature;  however,  at  11°C  the 
mirrors  remained  free  of  condensation.  Sons  of  the  data 
scatter  was  caused  by  drift  effects  of  the  spectroneter . 

The  calibration  of  a;  with  AT  (see  Figure  4)  determines  the 
dR- scale  and  la  consisted  by  a  measurement  of  the  Q  value 
via  the  bandwidth  of  the  resonance. 


Results  and  Discussion 

Experimental  Results 

Pressure  scans  of  the  attenuation  rate  n*  due  to  H?0 
absorption  were  made  at  the  frequency  ft  *  138.2  CHx. 

First,  pure  water  vapor  pressure  e  wee  increased  to  a 
constant  relative  tumidity  RH,  adjusted  for  values  between 
80  end  100  percent.  Then  nitrogen  wee  Injected  to  simulate 
■Diet  air,  although  Hj  broadens  water  vapor  lines  1.096 
•ore  effectively  then  dry  air  (15).  The  teat  frequency  f: 
falla  into  an  EKF  window  where  little  attenuation  by  local 
lines  le  expected  (aee  Eq.27).  A  sumnery  of  reduced  end 
calibrated  date  is  shown  In  Figure  6.  A  data  fit  to 
Eq. (13)  at  300  K  amounted  to 

J  •  0.197**  *  0.00750eP  dB/kn.  ':i) 


resulting  in 


26.3  •  f.  5n 


The  broadening  etf iciencv  Is  determined  by  pressure 
pairs  e:  and  points  winre  water  vapor  attenuation 
(*  c^)  and  moist  nitrosun  attenuation  (e  P(e-l)]  are  equal, 

»*  -  (F|/ei)  -1 

Attenuation  at  points  A  (2.60  d3/kn ,  e  •  3.49  k»a)  and  B 
(5.10  dB/ha,  P  •  106.7  kPa'  was  studied  with  special  care. 
Point  A  was  stable  over  a  -0-honr  period  and  independant 
of  the  water  valve  position  lopen/closed),  provided  a 
•light  mirror  heating  (IV)  was  applied.  Increasing  the 
mirror  heating  to  5W  (+4  C)  did  not  change  the  result. 

Strong  anomalies  (condensation  effects),  on  the  other 
hand,  were  observed  when  the  mirror  heat  era  were  off.  At 
point  8  e  very  long  mixing  tire  (several  hours)  lies  to  be 
allowed  when  the  water  valve  le  opened  against  110  kPa 
nitrogen.  Fester  nixing  occurs  when  nitrogen  is  blended 
with  water  vapor  alreedv  filling  the  cell. 

ijesl  Tie*  S|«t  ♦  *tr  7  t  *••!  »  i  ••  •!<**  e»r  7  1MJ 

PICS  ZOO  T2_  ??  .  T*  J  .  S  J  3  M  M  f4  >1  > 


1SUN  S  1  TOPE  S  20 


son  t  «  raoi  e  *e 


•  !<•)  «  fk  0.3041/h.  y«un.<l  y  = 

l  •  •  1.102  kPa  \  1  °t.- 


i»  A. > . 1 1  isfj'iivTi 
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Figure  5.  Example  of  a  ties  series  date  plot  froa  five 
lnscraeants,  reeding  temperature:  Tj  •  alrror,  T»,a  - 
cell  center  end  bottom,  Qj  •  water;  end  peak  signal  Pk 
(water  vapor  ebeorptlan).  The  sequence  covers  e  two-hour 
period  during  which  wxter  vapor  saturation  was  achieved  by 
equalising  the  teaperatura  of  8.-  C  between  cell  (T*,«) 
i.nd  water  (<h).  CAL  -  calibration  points  for  *\  (dB),  BV 
(■•)  *  reference  power,  m4  nl  crows  tar  setting. 


Temperature  wupenuenre  la  another  diagnostic  which 
can  bu  used  to  support  tn*  theoretical  uuerpretat ion  of 
continuum  attenuation  in  tiros  oi  nolcrular  line 
broadening,  a  series  of  pure  water  vapor  pressure  scans 
was  made  at  T  •  282  K  (see  figura  5).  Temperature  limited 
the  saturation  pressure  to  1.132  i2a.  A  preliminary  fit 
of  these  data  yields  (coapare  with  Eq.15) 


Pressure  dependence  indicated  by  Eq.(13)  has  been  ob¬ 
served  in  earlier  laboratory  (10  to  17,  23,  24)  and  field 
(4d,  4e)  experiments.  A  selected  review  of  laboratory 
data  is  given  in  Table  IV.  The  following  approximate  re¬ 
latione  can  be  deduced  from  this  decs  body: 

a,  *  hgte2  ♦  ep/m*)  (la  between  llaes); 

kft  «  f :  (excluding  the  67.7  CHx  data); 
e*  •  l/f°-S  (m.  Uf.  10); 

j  ’  5.5  (for  th.  110  GKz  dec*). 

Absolute  date  on  are  difficult  to  obtain.  In  spite  of 
existing  inconsistencies,  our  data  lend  support  to  the 
hypothesis  that  far-wings  of  self -broadened  H*0  lines 
centered  above  1  THx  ere  probably  the  sajor  source  for  ehe 
EHF  continuum  (10).  Line  *  hap  e_  theory  neads  to  be  revised 
to  eeeount  for  the  minute  (<  10‘-  the  value  at  line  center) 
contributions  in  the  extreme  low  frequency  wings  of 
prsssure-broadened  lines  located  above  l  THx  in  the  rota¬ 
tional  H:0  spectra  (4b,  4h,  9,  21,  22). 

Table  IV.  Laboratory  results  for  water  vapor  attenuation 
a  fitted  to  Eq.(13)  yielding  the  coefficient  k*  end 
broadening  efficiency  nx.  s 


f 

k* 

a 

X 

n 

Broadening 

Atent 

T 

Ref. 

CHx 

dB/Ckm  kPa 

*> 

t 

31 

0.011 

49* 

Air 

318 

U 

67.7 

0.1* 

44* 

Air 

283 

24 

110 

0.12 

35* 

Air 

304 

2* 

UO 

0.16 

35* 

Air 

292 

23 

110 

0.21 

35* 

Air 

274 

2* 

13B 

0.197 

26.3(29) 

::,fAir) 

300 

this 

138 

0.25 

282 

paper 

213 

0.62 

24* 

?;2 

320 

4g,  17 

3*5 

2.4  :  1.5 

24 

300 

19 

ill 

2.9  :  0.8 

18 

“2 

300 

19 

667 

4.1  :  1.1 

9.6 

*•** 

300 

19 

883 

5.2  1.4 

9.4  ±  0. 

.9 

300 

19 

Our  fit  to  reported  raw  data 


Moist  Micro  gen 
136.2  CHx 
300  K 
(•#-3.49) 

282  K 

(«s«1.10) 


0  12  3  4  <•» 

Vapor  Fra.sur*  •  jg  73  100 

TOTAL  PRESSURE  F,  kP* 

figure  6.  Yrtuurt  dependence  of  specific  ,tt«watloD 
a  duo  to  water  vopor  (o)  and  moist  nitrogen  (P)  aaa- 
aured  at  138.2  CHi  and  relative  humidities  fro*  *H  -  80 
to  100  percent. 
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Water  Vaoor  Continuum  and  E!?F  Preparation  ’todellrg 

The  experimental  result  at  138.2  CHz  for  colst  nitro¬ 
gen  attenuation  (Eq.24)  has  been  incorporated  into  the  EHF 
propagation  model  chat  ia  baaed  upon  LB 3  (1,23).  Reducing 
che  ep-ccrm  by  9.4  percent  for  air-broadening  and  aubtraet- 
ing  line  contribution*  due  to  LB5, 

12(US)  -  l.f*  •  10'3  f  *  0.0093  «*  dl/kB.  (27) 

gives  the  pressure  dependence  of  continuum  attenuation  as 

-  *.»6  •  10*3  tp  ♦  0.ua«:  dS/ta.  (28) 

Adding  frequency  and  reaper ature  dependences  to  Eq .29 
leads  to  a  moist  air  continuum  spectrum  for  LS5t  modeled  by 

!T  -  (cl«p6)'  ♦  c,«V*U  »  10'6  ppo.  <291 


The  coefficients  for  Eq.(29)  are  chosen  according  to  the 
following  achm: 


C1 

y 

c3 

y* 

Remark 

Old 

1.90 

3.1 

0 

l 

vUuc-Reifensceia 

Vcw 

1.40 

2.5 

54.1 

3.3 

Fes. 28.  13.  18 

The  "new”  continuum  spectrum  is  different  fron  various  other 
empirical  forms  (summarized  in  Table  V)  which  have  been  pro¬ 
posed  . 

Table  VI  compares  che  implications  of  the  old  and  the 
new  continuum  for  the  case  of  predicting  xenlta  attenua¬ 
tion  MdB)  through  a  model  atmosphere. 


Absolute  attenuation  rate*  at.  have  been  -assured  in 
the  laboratory  for  pure  water  vapor  and  tor  fixtures  of 
vapor  end  nitrogen,  sinw-atir.t  moist  air.  Initial  studies 
were  performed  at  128  G’iz,  two  temperatures  (292  anc 
300  K).  partial  vapor  pressures  up  to  saturation  (1.10  and 
3. 48  kPa),  and  total  pressures  up  to  160  kPa.  Preliminary 
results  were  acquired  for  pressure  and  temperature  depend¬ 
ences  of  ay.  The  data  behaved  ia  a  fashion  similar  to 
that  expected  from  H?0  lias- broadening  theory  (Eq.L3). 
Except  for  a  strong  self -broadening  component  («  e~), 
no  additional  anomalous  behavior  was  observed— even  at 
saturation  pressures  (RB  •  100S) — contrary  to  findings  by 
other  workers  (4,17,  etc.). 

Zn  our  case,  wa  could  clearly  identify  ssnmalnue  ab¬ 
sorption  behavior  (a.g. ,  unusually  high  rates  r  extreme 
toftmure  dependences,  hyataraaea  ia  pressure  and  tam- 
paratura  cycles)  as  being  Instrumental,  that  is  condensa¬ 
tion  effects  on  the  millimeter  wave-active  parts  of  che 
spectrometer.  Considerable  effort  had  to  be  expended 
before  the  inscnoMtnt  produced  consistent,  reproducible 
reeults . 

A  water  vapor  cencisutmi  N”  was  formulated  for  the  EHF 
propagation  model  (1,25)  in  order  to  supplement  contribu¬ 
tions  from  30  B?0  lines  below  1  THz-  At  138  CHz ,  the 
local  lines  furnish  40  percent  to  eir- broadening  but  only 
3  percent  to  self- broadening  cf  the  water  vapor  attenua¬ 
tion  Oy.  Practical  laplicatlcr.s  of  Che  continuum  to 
propagation  modal  mg  are  not  t'a  serious.  Up  to  humidi¬ 
ties  of  v  4  10  g/s^,  the  widely  used  empirical  Caut- 
Relfenstein  continuum  <"1/  is  roughly  confined.  At 
higher  values.  Increases  over  che  CH  continuum  and  a 
nonlinear  dependence  on  absolute  faeldlty  sre  predicted. 
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table  V.  Experimental  water  vapor  continuum  coefficients 
(Eq.6)  reported  by  other  workers  (2). 


k 

«<") 

t(T) 

o(» 

Local  Linas 

f.  (Ex 

Raf . 

0.023 

1.00 

2.1 

2.50 

HVW 

50  to  500 

10 

0.039 

1.05 

3.0 

2.40 

GR 

20  to  350 

18 

0.049 

1.00 

2.1 

2.00 

vw 

15  to  500 

<*) 

0.061 

1.00 

2.1 

1.22 

gr 

150  to  410 

21 
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Table  VI.  Comparison  of  zenith  attenuation  A(dB)  predicted 
for  a  synthetic  atmosphere  (30°H,  July)  by  employing  the 
E8F  Propagation  Model  (1, 25)  with  old  (Cant-tslfanstein) 
and  new  (Eq.29)  water  vapor  continuum  formulations. 
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